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NOTICES 


Disclaimers 

The  findings  in  this  report  are  not  to  be  construed  as  an 
official  Department  of  the  Army  position,  unless  so  desig* 
nated  by  other  authorized  documents. 

The  citation  of  trade  names  and  names  of  manufacturers  in 
this  report  is  not  to  be  construed  as  official  Government 
indorsement  or  approval  of  commercial  products  or  services 
referenced  herein. 
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INTRODUCTION 

The  purpose  of  smoke  screening/obscuration  with  chemical  smokes  such  as 
white  phosphorus  (WP)  or  hexachloroethane  (HC)  in  a tactical  environment 
is  to  attenuate  the  transmission  of  electromagnetic  radiation  in  a finite 
optical  path  or  line  of  sight  to  a threshold  level  of  the  human  eye  or 
electro-optical  sighting,  designating,  or  guidance  system.  Attenuation 
of  electromagnetic  radiation  is  dependent  upon  such  factors  as  absorption 
and  scattering  i9««MaaaaiflHMMl[BiMcaipi t»tion  and  the  crosswind  integrated 
concentration  of  the  plume  or  puff  of  smoke  deliberately  injected  into  the 
line  of  sight.  The  amount  of  chemically  produced  smoke  (munition  expen- 
ditures) required  can  be  determined  from  empirically  derived  relationships 
of  transmittance  versus  concentration  data. 

Current  techniques  for  calculating  munition  expenditures  do  not  consider 
certain  meteorological  parameters  such  as  visibility  and  relative  humid- 
ity. Ignoring  such  parameters  can  result  in  the  incorrect  calculation  of 
required  smoke  munitions.  Consequently,  a semi  empirical  algorithm  has 
been  developed  for  battlefield  use  for  calculating  munition  expenditures. 

The  algorithm  is  called  KWIK  (an  acronym  for  crosswind  integrated  concen- 
tration). The  primary  output  is  munition  expenditures  for  artillery  en- 
gaged in  smoke  screening/obscuration  support  of  other  military  elements. 
This  report  describes  the  meteorological  and  mathematical  background  of 
KWIK. 

Figure  1 shows  a typical  scenario  where  smoke  screening  is  utilized  to 
conceal  the  movement  of  friendly  forces  along  path  X.  The  enemy  is  be- 
lieved to  be  somewhere  on  a hill  having  a line  of  sight  L.  The  symbols 
used  on  the  figure  are  described  as  follows: 

N - grid  north 

S - angle  of  line  of  sight  to  target 
L - slant  range  to  target 
X - distance  to  be  smoked 
Z - mean  height  of  target 
H - release  height  of  smoke 
V - wind  direction 

oc  - direction  of  line  of  sight  from  grid  north 

A - angle  between  wind  direction  and  line  of  sight:  |V  - a| 

A block  diagram  is  presented  in  figure  2 which  contains  details  of  the 
four  main  sections  of  the  KWIK  smoke  program:  (1)  meteorological  inputs 
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and  calculations,  (2)  meteorological  optics  and  smoke  concentration  cal- 
culations, (3)  atmospheric  diffusion  and  smoke  source  strength  calcula- 
tions, and  (4)  munition  expenditures.  A description  of  the  algorithm  is 
presented  to  show  mathematically  how  the  four  parts  of  the  algorithm  are 
derived.  The  fiftegi  tabies  contain  equations  and  values lOf  coefficients 
which  are  essential  to  the  appropriate  calculations  in  the  algorithm.  A 
flowchart  of  the  algorithm  is  also  given  followed  by  glossary  of  mnemonics 
aiwK  software  ttstings  u#t7ig  FORTRAN,  BASrC,  AND  HPL  programming  lalf^uages. 

» The  latter  is  for  use  with  the  Hewlett-Packard  model  9825A  desktop  computer. 


DESCRIPTION  OF  ALGORITHM 

Meteorological  Inputs  and  Calculations 

Nine  meteorological  inputs  are  necessary  to  execute  the  algorithm.  The 
inputs  are: 

1.  Cloud  ceiling  (hundreds  of  feet) 

2.  Cloud  cover  (percent) 

3.  Visibility  (miles) 

4.  Precipitation  indicator  (yes  or  no) 

5.  Temperature  (degrees  F) 

6.  Dewpoint  (degrees  F) 

7.  Wind  direction  (tens  of  degrees) 

8.  Windspeed  (knots) 

9.  Length  of  average  surface  roughness  element  (centimeters) 

Met  information  of  this  type  is  available  in  the  battlefield  from  one  of 
the  following  sources: 

1.  US  Army  observation 

2.  US  Air  Force  Air  Weather  Service  (USAF-AWS)  airfield  observations 

3.  AF-AWS-GWC  (Global  Weather  Central)  prognostications 

In  addition,  the  following  information  about  the  location  requiring  the 
smoke  is  also  an  input: 

1.  Site  identification  (if  any)  i 

2.  Latitude  of  site  (decimal  degrees) 


3.  Longitude  of  site  (decimal  degrees) 

4.  Altitude  of  site  (kilometers) 

5.  Julian  data  and  Zulu  hour  of  met  data  recording 

6.  Slant  range  to  target  (meters) 

7.  Angle  of  sight  to  target  (decimal  degrees) 

8.  Total  distance  to  be  smoked  (meters)  * 

9.  Release  height  of  smoke  (meters) 

10.  Mean  height  of  target  (meters) 

11.  Direction  of  line  of  sight  (decimal  degrees) 

12.  Time  smoke  is  required  (minutes) 

The  Input  units  are  shown  both  in  the  English  and  metric  systems  because 
that  is  the  form  in  which  ttiey  are  supplied  by  the  sources  quoted.  However, 
the  program  converts  all  units  to  the  metric  system. 

The  calculation  of  atmospheric  stability  is  based  upon  work  by  Dr.  F. 
Pasquillo^  Stability  close  to  the  ground  is  dependent  mainly  upon  net 
radiation  and  windspeed.  Incoming  radiation  is  dependent  upon  solar 
altitude,  which  is  a function  of  time  of  day  and  season  of  the  year. 

The  amount  of  cloud  cover  and  its  thickness  will  also  influence  incoming 
or  outgoing  radiation.  For  daytime,  table  1 is  used  to  arrive  at  an 
insolation  class  number  as  a function  of  solar  altitude.  This  number 
becomes  the  net  radiation  index  (table  2)  after  being  modified  by  the 
amount  of  cloud  cover  and  ceiling.  For  instance,  if  the  total  cloud  cover 
is  100  percent  and  the  ceiling  is  less  than  7000  feet  (whether  day  or 
night),  the  net  radiation  index  is  equal  to  0.  For  nightime,  estimates 
of  outgoing  radiation  are  made  by  considering  the  amount  of  cloud  cover. 

For  example,  if  the  total  cloud  cover  is  less  than  or  equal  to  40  percent, 

-2  will  be  used  for  the  net  radiation  index;  if  it  is  greater  than  40 
percent,  the  net  radiation  index  will  be  -1.  Table  2 shows  the  stability 
class  as  a function  of  windspeed  and  radiation  index.  The  stability 
classes  are  identified  as  follows: 

A - extremely  unstable 

B - unstable 

C - slightly  unstable 


^F.  Pasquill,  1961,  "The  Estimation  of  the  Dispersion  of  Windborne  Material," 
Meteorol  Mag,  Vol  90 
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0 - neutral 
E - slightly  stable 
F - stable 

G - extremely  stable 

These  are  also  identified,  for  computation  purposes,  as  numbers  1 through 
7.  The  net  radiation  index  ranges  from  4,  the  highest  positive  net  radia- 
tion directed  toward  the  ground,  to  -2,  the  highest  negative  net  radiation 
directed  away  from  the  earth.  Instability  occurs  when  the  positive  net 
radiation  is  high  and  winds  are  light  and  during  neutral  conditions  with 
cloudy  skies  or  high  windspeeds.  The  algorithm  utilized  this  procedure 
to  determine  the  atmospheric  stability  category. 

The  algorithm  converts  visibility  from  miles  to  kilometers  and  temperature 
and  dewpoint  from  degrees  Fahrenheit  to  degrees  Celcius.  Relative  humidity, 
a function  of  temperature  and  dewpoint,  is  calculated  as  a ratio  of  the 
vapor  pressure  (P^)  at  a temperature  T (degrees  C)  and  the  saturation  vapor 

pressure  (Pyj^)  at  dewpoint  T^  (degrees  C)  expressed  by: 

P 

RH  = -^  (100)  , (1) 

P 

vm 


I 


where 


vm 


= 6.11(10)  ° ° 


The  constants  a and  b are  as  follows: 
over  ice 
a = 9.5 
b = 265.5 


over  water 
a = 7.5 
b = 237.3 


The  above  solution,  according  to  Haurwitz^  can  be  attributed  to  Tetons  and 
is  similar  to  a theoretical  formula  derived  from  the  equation  of  Clausius 

1 


^B.  Haurwitz,  1941,  Dynamic  Meteorology,  McGraw-Hill  Book  Company,  Inc., 
New  York 
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Clapeyron  described  by  Brunt. ^ The  only  condition  necessary  to  apply 
the  above  is  the  ability  to  distinguish  between  the  saturation  pressure 
over  ice  and  over  water,  i.e.,  the  freezing  point. 


Meteorological  Optics  and  Smoke  Concentration  Calculations 

The  optics  portion  of  KWIK  is  adapted  from  an  approach  to  atmospheric 
transmission  suggestecTby  Downs.**  Transmittance  of  light  at  various 
wavelengths  through  a path  is  determined  by  calculating  the  attenuation 
due  to  absorption  by  water  vapor,  scattering  by  haze  or  fog,  and  precip- 
itation. When  the  attenuation  due  to  atmospheric  conditions  is  known, 
the  attenuation  due  to  smoke  that  is  required  to  lower  transmittance  to 
a threshold  contrast  for  a particular  wavelength  can  be  computed.  By 
use  of  the  transmittance  and  empirically  derived  relationships  between 
transmittance  and  concentration  for  various  smokes,  the  crosswind  inte- 
grated concentration  for  a particular  smoke  can  then  be  computed. 

Absorption  is  directly  attributable  to  the  amount  of  preci pi  table  water 
in  a path,  assuming  the  water  vapor  concentration  in  the  atmosphere  is 
reasonably  well  behaved  and  exhibits  a scale  height  of  about  2 km.  The 
water  vapor  concentration  expressed  in  centimeters  per  kilometer  of  path 
length  may  then  be  given  as: 


0 


(2) 


( 


\ 


where  W^  is  the  precipitable  water  along  a path  L and  0 is  the  angle 

between  the  horizontal  and  the  height  of  a target  above  or  below  an  obser- 
ver. is  computed  from  the  following  linear  regression  equation  re- 
lating precipitable  water  and  dewpoint  temperature  (T^): 


W^  = 0.4477  + 0.0328T^  + 1.2(10)"3t2  + 1.84(10)"5t3  (3) 

Equation  (3)  was  fit  to  data  extracted  from  Downs'*  (figure  4 in  Downs)  , , 

and  is  considered  valid  for  all  geographical  regions.  ! 


3d.  Brunt,  1952,  Physical  and  Dynamical  Meteorology,  second  edition,  Cambridge 
University  Press,  London 

I ‘*A.  R.  Downs,  1976,  "A  Review  of  Atmospheric  Transmission  Information  in 

i the  Optical  and  Microwave  Spectral  Regions,"  Ballistics  Research  Laboratories 

I Report  2710 
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The  amount  of  water  vapor  in  the  path,  W,  is  given  by: 

W - 


W - Wq  /*  dL  . 

» *'o»  » 


(4) 


Transmission  through  the  absorbing  component  of  the  atmosphere  is  calcu- 
■'MALUU  by  using  an  error  function  Absorption  law  developed  by  Elsasser^ 


T = 1 - erf  (z)  , 

(5) 

^2 

(6) 

where 


z = 0. 5 3 /ttW 


3 = error  function  absorption  coefficient  as  a function  of  wave- 
length (table  3). 

Downs'*  states  that  the  Elsasser  approach  is  unable  to  correctly  address 
long  wavelengths  and  suggests  using  an  approach  described  by  Fisher®  for 
the  far  infrared  wavelengths  (8vm-14ym).  Thus  the  computation  of  trans- 
mission due  to  absorption  by  water  vapor  for  the  long  wavelength  case  is 
given  by: 


j , g-0.0681W 


(7) 


Reduction  in  transmittance  due  to  attenuation  by  haze  and  fog  can  be 
calculated  by  using  the  Mie  theory.  Downs'*  indicates  that  the  Mie 
scattering  coefficient  decreases  with  altitude  such  that  its  behavior 


®W.  M.  Elsasser,  1942,  "Heat  Transfer  by  Infrared  Radiation  in  the 
Atmosphere,"  Harvard  Meteorological  Series  6,  Harvard  University  Press, 
Cambridge,  MA 

'*A.  R.  Downs,  1976,  "A  Review  of  Atmospheric  Transmission  Information 
in  the  Optical  and  Microwave  Spectral  Regions,"  Ballistics  Research 
Laboratories  Report  2710 

®D.  F.  Fisher  et  al . , 1963,  "Transmissometry  and  Atmospheric  Transmission 
Studies  Final  Report,"  University  of  Michigan,  Institute  of  Science  and 
Technology 
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can  only  be  estimated.  Thus  the  following  expressions  for  (Mie 

scattering  coefficient)  are,  at  best,  an  approximation  to  the  behavior 
of  the  Ou  versus  altitude  relationship 

„ „-L  sin  0/4.1  M ^ iQ\ 

= Owf®  V > G{X)  , (8) 


" ^hf®*"  ° y ^ g 0 < L sin  0 < 1 km 

= 0.128e"'-  V - G (A),  1 km  < L sin  0 < “ 

M2 


o 


M 


(9) 


where  is  a coefficient  of  attenuation  determined  from  a linear  regression 

equation  as  a function  of  visibility  and  wavelength  (table  4),  based  upon 
Downs'**  evaluation  of  the  best 'available  data.  V is  visibility,  and  G(a) 
is  the  scale  height  of  G(X)  is  not  constant;  rather  it  is  a function 

of  altitude,  visibility,  and  wavelength.  Table  5 indicates  approximations 
for  G(A)  for  various  wavelengths  considered.  Transmission  along  a path 
with  attenuation  can  be  determined  by  the  equation 


T - 


(10) 


Thus,  reduction  in  transmittance  due  to  attenuation  by  haze  and  fog  can 
be  calculated  by  using  equation  (10)  by  substituting  a value  for  accord- 
ing to  equations  (8)  and  (9).  If  precipitation  is  Indicated  (by  input 
parameter),  then  the  value  for  transmittance  in  equation  (10)  is  set  to 
one  and  a calculation  is  made  for  attenuation  by  precipitation  instead. 


“♦A.  R.  Downs,  1976,  "A  Review  of  Atmospheric  Transmission  Information 
in  the  Optical  and  Microwave  Spectral  Regions Ballistics  Research 
Laboratories  Report  2710 
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The  reduction  in  transmittance  due  to  attenuation  by  precipitation  can  be 
obtained  from  the  equation 


where  L is  the  path  length  and  is  an  attenuation  coefficient  deter- 
mined from  a linear  regression  equation  as  a function  of  visibility  and 
wavelength  (table  6). 

The  total  transmittance  along  an  optical  path  is  the  product  of  the 
partial  transmittance 


^total  ~ ^a^hf^p^s 


(12) 


where 


T = Transmittance  due  to  attenuation  by  atmospheric  absorption 

d 

T^^  = Transmittance  due  to  attenuation  by  haze  and  fog 

Tp  = Transmittance  due  to  attenuation  by  precipitation 

T = Transmittance  due  to  attenuation  by  smoke 
s 

T can  then  be  calculated  from  equation  (12),  and  the  desired  thresh- 
s 

hold  contrast  of  transmittance  for  a particular  wavelength  can  be  expressed 
by  using  the  following  equation: 


T 

s 


^a^hf'’’p 


(13) 


where  T^^,  the  threshold  contrast,  is  based  upon  the  Koschmieder'^  theory 
and  is  set  equal  to  0.02. 

After  the  transmittance  due  to  attenuation  by  smoke  has  been  computed,  the 
line  of  sight  integrated  concentration  (CL)  necessary  to  achieve  this 
value  can  be  calculated  from  a linear  regression  equation  as  a function 
of  transmittance  and  wavelength  (table  7).  The  data  used  to  fit  the 
regression  equations  comes  from  laboratory  and  atmospheric  test  results 
for  transmittance  through  variable  concentrations  of  different  smokes 
over  finite  path  lengths.  The  line  of  sight  concentrations  may  now  be 
utilized  in  a Gaussian  diffusion  scheme  to  determine  smoke  source  strengths. 


^H.  Koschmieder,  1924,  "Theorie  der  Horizontalen  Sichtweite,"  Beitr 
Phys  Frein  Atmos,  12:33-53,  171-181 
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Batchelor®  pointed  out  that  the  Gaussian  function  could  provide  a general 
description  of  average  diffusion  in  a continuous  plume.  Diffusion  studies 
by  Hay  and  Pasquill;®  Cramer,  Record,  and  Vaughan;^®  and  Barad  and  Haugen^^ 
suggested  that  Gaussian  plume  formulae  are  quite  practical  and  applicable 
in  the  atmosphere. 


Atmospheric  Diffusion  and  Smoke  Source  Strength  Calculations 

Successful  smoke  screening/obscuration  in  the  surface  boundary  layer  is 
largely  dependent  upon  the  existence  of  the  meteorological  data  necessary 
to  describe  the  diffusion  process.  Generally,  because  of  the  type  of 
meteorological  observations  available  in  a battlefield  environment,  an 
assumption  of  Gaussian  diffusion  may  be  the  only  model  for  which  adequate 
data  will  exist.  Since  KWIK  was  developed  for  actual  field  use,  Gaussian 
diffusion  theory  was  employed.  Estimates  of  downwind  dispersion  may  be 
derived  by  using  methodology  developed  by  various  individuals,  including 
F.  Pasquill,  F.  A.  Gifford,  F.  B.  Smith,  D.  B.  Turner,  and  R.  P.  Hosker,  Jr. 

The  vertical  dispersion  coefficient,  o , is  calculated  from  the  follow- 
ing equation:  ^ 


(14) 


> 10  cm 


< 10  cm 


o,  = F{2^;  X)  g(X) 

U 


'"i 

|c.  X 

|1  + (C2  X 

F(2  ; X)  = 

L J 

In , 
1 

[{1  + C2  X 

and 


g(X)  = (ai  X ^i)/(l  + aa  X ^2)  , 


®G.  K.  Batchelor,  1949,  "Diffusion  in  a Field  of  Homogeneous  Turbulence, 

I.  Eulerian  Analysis,"  Australian  J Sci  Res,  2:437-450 

®J.  S.  Hay  and  F.  Pasquill,  1957,  "Diffusion  from  a Fixed  Source  at  a 
Height  of  a Few  Hundred  Feet  in  the  Atmosphere,"  J Fluid  Mech,  2:299 

i°H.  E.  Cramer,  F.  A.  Record,  and  H.  C.  Vaughan,  1958,  The  Study  of  the 
Diffusion  of  Gases  or  Aerosols  in  the  Lower  Atmosphere.  Report  AFTRC-TR- 
58-239,  Department  of  Meteorology,  Massachusetts  Institute  of  Technology 

^^M,  L.  Barad  and  D.  A.  Haugen,  1959,  "A  Preliminary  Evaluation  of  Sutton's 
Hypothesis  for  Diffusion  from  a Continuous  Point  Source,"  J Meteorol , 16:12 
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where 

The  coefficients  ai,  bj,  a2,  b2  are  given  in  table  8, 
the  coefficients  Cj,  di,  C2,  d2  are  given  in  table  9, 
and 

X is  total  distance  to  be  smoked  (meters). 

is  the  surface  roughness  length  computed  by  using  the  average  rough- 
ness element  (Y  - one  of  the  program  inputs)  and  the  following  equation: 


logio  ^ -1.24  + 1.19  log^o  Y (15) 

where  Y is  estimated  by  a visual  survey  of  the  average  height  of  trees, 
bushes,  or  grass. 

The  lateral  dispersion  coefficient,  a^,  is  calculated  from  the  follow- 
ing equation:  3 


ay(X)  = a X^  (16) 

where  X is  the  total  distance  to  be  smoked,  a is  obtained  from  table  10 
as  a function  of  stability  category,  and  b is  a constant  equal  to  9.0. 

The  previous  calculations  of  and  are  for  a continuous  source,  such 

as  HC  smoke.  For  a quasi-instantaneous  source,  such  as  WP,  and 

are  computed  as  two- thirds  of  the  values  determined  in  the  continuous  case, 
based  on  experimental  results  of  Smith  and  Hay.^^ 

The  crosswind  integrated  concentration,  CWIC,  is  determined  from  the  computed 
line  of  sight  integrated  concentration  value  by  applying  a wind  direction 
correction  factor,  f,  as  follows: 


CWIC  = (CL)f  , 


(17) 


B.  Smith  and  J.  S.  Hay,  1961,  "The  Expansion  of  Clusters  of  Particles 
in  the  Atmosphere,"  Quart  .1  Roy  Meteorol  Soc,  87:82 
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where 


f = /3. 31 24/ [6. 76  (sin  A)2  + 0.49  (cos  A)2J 
CL  = line  of  sight  integrated  concentration 
A = Ml  - d2| 
di  = surface  wind  direction 
d2  = direction  of  line  of  sight 


After  the  dispersion  coefficients  and  the  crosswind  integrated  concentra- 
tion have  been  determined,  source  strength,  Q,  can  be  computed  from  the 
following  equations: 

Continuous  Source 


(18) 


Quasi -Instantaneous  Source 


Q 


(19) 


where 


0^  = vertical  dispersion  coefficient 
= lateral  dispersion  coefficient 
H = release  height  of  smoke  source  (meters)  (figure  1) 

Z = mean  height  of  target  (meters) 

V = surface  windspeed  (mps) 

Chemical  smoke  aerosols  are  hygroscopic  and  thus  the  source  strength  must 
be  modified  by  a yield  factor,  which  is  a function  of  relative  humidity, 
as  follows: 


= Q/Y  . (20) 

where  Y is  the  yield  factor  (table  11),  and  Op,  is  the  modified  soufce 
strength. 


Munition  Expenditures 


KWIK  calculates  munition  expenditures  following  the  procedure  described 
in  Army  Training  Circular  6-20-5. Basically,  this  procedure  determines 
the  number  of  guns  and  the  number  of  munitions  per  gun  in  order  to  calcu- 
late the  total  number  of  munitions  needed  for  a smoke  mission. 

The  number  of  guns  is  determined  by  the  following  equation: 


G = Q^/S 


(21) 


where 


G = number  of  guns 

S = unit  source  strength  (i.e.,  total  amount  of  chemical  material 
available  in  one  munition) 


The  number  of  guns  calculated  is  always  rounded  up  to  the  nearest  whole 
number. 

Next,  the  total  time  for  replenishment  (i.e.,  the  total  time  smoke  replen- 
ishment is  required  to  maintain  the  desired  screen)  is  obtained  from  the 
relation: 


W = T + A - B (22) 

where 

W = time  for  replenishment 

T = total  time  smoke  is  needed 

A = time  required  for  buildup  (table  13) 

B = munition  burn  time  (table  14) 

Since  time  for  replenishment  has  been  calculated,  the  number  of  munitions 
required  for  one  gun  can  be  computed  as  follows: 


F = 1 + WD 


(23) 


^^"Field  Artillery  Smoke,"  1975,  Training  Circular  6-20-5.  US  Army 
Field  Artillery  School,  Fort  Sill,  OK 


w 


where 


F = number  of  munitions  for  one  gun 
W = computed  time  for  replenishment 
D = rate  of  fire  of  gun  (table  15) 

The  constant  (1)  in  equation  (23)  represents  the  initial  munition  round 
fired  by  a gun. 

Total  munitions  required,  R,  can  now  be  calculated  by  the  following 
aquation: 


R = GF 


(24) 


where 


G is  the  number  of  guns,  and  F is  the  number  of  munitions  per  gun. 


M 


Smoke  sereenini^  scenario 


KWIK  SMOKE  PROGRAM 


TABLE  1.  INSOLATION  AS  A FUNCTION  OF  SOLAR  ALTITUDE 


Solar  Altitude 


Insolation 


cc 

Insolation 

Class 

o 

o 

< OC 

Strong 

4 

35° 

< a < 60° 

Moderate 

3 

15° 

< OC  < 35° 

Slight 

2 

OC  < 15° 

Weak 

1 

TABLE  2.  STABILITY  CLASS  AS  A FUNCTION  OF  NET  RADIATION  AND  WINDSPEED 


Winds peed 
(knots) 


Net  Radiation  Index 


-1  -2 


B B 
B B 


D D 

D D 


TABLE  3.  ERROR  FUNCTION  ABSORPTION  COEFFICIENTS 


Wavelength 

8 

Visual 

0.118 

1 . 06ym 

0.22 

2.3ym 

0.14 

3.8ym 

0.55 

TABLE  4.  HAZE  AND  FOG  ATTENUATION  COEFFICIENT  EQUATIONS 

Visual  In  = 1.5551  - 0.9811  InV  - 0.0197  {lnV)2  + 0.0041  (lnV)3 

1.06ym  In  = 1.5551  - 0.9811  InV  - 0.0197  (lnV)2  + 0.0041  (lnV)3 

nf 

2.3ym  In  = 1.4491  - 1.0044  InV  - 0.12(lnV)2  + 0.0032  (lnV)3 

hr 

3.8vitn  In  = 1.2394  - 1.0436  InV  + 0.0099  (lnV)2  - 0.0016  (lnV)3 
hf 

1.06yni  In  = 1.5176  - 1.7147  InV  + 0.001  (lnV)2  + 0.0428  (lnV)3 
ht 

= Haze  and  fog  attenuation  coefficient 
V = Visibility  (km) 


I 


TABLE  5.  G(X)  - APPROXIMATIONS  FOR  SCALE  HEIGHT  OF  HAZE  AND  FOG 
ATTENUATION  COEFFICIENTS  AS  A FUNCTION  OF  WAVELENGTH 


Wavelength 

G(X) 

Visual 

26.7 

1 .06vim 

13.0 

2.3ym 

5.3 

3.8ym 

5.1 

10.6yni 

5.0 

TABLE  6.  RAIN  ATTENUATION  COEFFICIENT  EQUATIONS 

Visual  In-  a = 1.3306  - 0.8825  InV  - 0.0753  (InV)^  + 0.0129  (lnV)3 
r 

l,06ym  In-a^  = 1.4098  - 0.9865  InV  - 0.0140  (InV)^  + 0.0023  (lnV)3 

2.3ytn  In- a = 1.5497  - 0.8696  InV  - 0.1084  (lnV)2  + 0.0231  (lnV)3 

r 

3.8uin  In- a = 1 .5556  - 0.9013  InV  - 0.0773  (lnV)2  + 0.0173  (lnV)3 
r 

10.6ym  In-a^  = 1.5928  - 0,9396  InV  - 0.0627  (InV)^  + 0.0168  (lnV)3 

= Rain  attenuation  coefficient 
r 

V = Visibility  (km) 


! 

i 
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TABLE  7.  EQUATIONS  FOR  DETERMINING  LINE  OF  SIGHT  INTEGRATED 

CONCENTRATION  (CL)  AS  A FUNCTION  OF  TRANSMITTANCE  (T) 
AND  TYPE  OF  CHEMICAL  SMOKE 


Visual 

Fog  oil  CL  = 0.0093  - 0.3428  InT  - 0.0009  (InT)^ 
HC  CL  = 0.0119  - 0.2747  InT  - 0.0013  (InT)^ 
FS  CL  = 0.0142  - 0.111  InT  + 0.00004  (lnT)2 
WP  CL  = 0.0055  - 0.1541  InT  - 0.0004  (InT)^ 


TABLE  8.  COEFFICIENTS  OF  THE  FUNCTION  g(X)  USED  IN  CALCULATING  THE 
VERTICAL  DISPERSION  COEFFICIENT  a2(X)  FOR  THE  VARIOUS 

STABILITY  CATEGORIES  (X  IS  GIVEN  IN  METERS) 


Stability  Category 

9l 

bi 

32 

t>2 

A 

0.112 

1.06 

5.38  (10-'+) 

0.815 

B 

0.130 

0.950 

6.52  (lO-**) 

0.750 

C 

0.112 

0.920 

9.05  (10-**) 

0.718 

D 

0.098 

0.889 

1.35  (10"3) 

0.688 

E 

0.0609 

0.895 

1.96  (10-3) 

0.684 

F 

0.0638 

0.783 

1.36  (10-3) 

0.672 

TABLE  9.  EQUATIONS  USED  TO  COMPUTE  THE  ROUGHNESS  CORRECTION 
FACTOR  F(z;  X)  USED  IN  CALCULATING 
a^CX)  [z^  IS  ROUGHNESS  ELEMENT  LENGTH  (cm)] 

In  Cl  = 0.444685869  + 0.294049265  (Inz^)  - 0.2  ’213914 
+ 0.155349504  (lnz^)3  - 0.032015723  (Inz^)** 

+ 2.15168  (10-3)  (lnz^)3 

In  di  = -1.298283909  - 1.006186784  (Inz  ) + 1.485094886  (Inz  )2 

- 0.774136725  (Inz  )3  + 0.156559355  (Inz  )*♦ 

0 0 

- 0.010823351  (Inz  )3 

0 

If  z^  < 10,  then 

In  C2  = 5.77267  (lO-**)  + 2.31943  (10-3)  (Inz^)  + 3.71041  (10-3)  (lnz^)2 

- 8.40602  (10-3)  (lnZjj)3  + 1.3421  (lO’^)  (Inz^)** 

+ 2.55131  (10-8)  (lnZjj)5 

If  10  £ Zp  £ 40,  then 

In  C2  = -11.56134901  = 2.148242814  (Inz^)  - 0.156210817  (lnZQ)2 
+ 7.03582  (10-3)  (inz^)3  - 1.47353  (10-‘*)  (lnz^)'+ 

+ 1.18256  (10-3)  (lnz^)5 

If  z >40,  then 
0 

In  C2  = 1108.366588  - 103.5495836  (Inz^)  + 2.424499256  (lnz^)2 

- 0.014584773  (Inz  )3  + 4.34517  (10-3)  (Inz  )*♦ 

0 0 

- 4.69556  (10-3)  (Inz  )3 

0 

In  d2  = 0.500775609  + 1.092615788  (Inz^)  - 1.573065836  (lnz^)2 

+ 0.724276579  (Inz  )3  - 0.140820904  (Inz  )*♦ 

0 0 

+ 9.61621  (10-3)  )5 

0 ' 


t 

I 


■ 


TABLE  10.  COEFFICIENT  USED  IN  CALCULATING  THE  LATERAL  DISPERSION 
COEFFICIENT  Oy(X)  FOR  THE  VARIOUS  STABILITY  CATEGORIES 

(X  IS  GIVEN  IN  METERS) 


Stability  Category  a 

A 0.22 
B 0.16 
C 0.11 
D 0.08 
E 0.06 
F 0.04 


TABLE  11.  EQUATIONS  USED  TO  CALCULATE  CHEMICAL  SMOKE  YIELD  FACTOR 
(Y)  AS  A FUNCTION  OF  RELATIVE  HUMIDITY  (RH) 

For  HC  Smoke; 

Y = 0.9337  + 0.0369  RH  - 7.0  (lO"*')  RH2  + 6.11  (10"^)  RH3 
For  FS  Smoke: 

Y = 1.3775  + 0.09868  RH  - 1.8  (10-3)  RH2  + 1.56  (10-3)  rh3 
For  Fog  Oil : 

Y = 1 

For  VIP  Smoke: 

Y = 3.2469  + 0.0774  RH  - 1 .6  (10-3)  rh2  + 1.73  (lO'S)  RH3 


TABLE  12,  UNIT  (PER  GUN)  SOURCE  STRENGTHS  (GRAMS) 


Gun 

105  Howitzer 

155 

Howitzer 

Munition 

HC  WP 

HC 

UP 

Source  Strength 

18.9  1737.3 

48.1 

5 7076.2 

TABLE  13. 

SMOKE  BUILDUP  TIME  (MIN) 

Gun 

105  Howitzer 

or 

155 

Howi tzer 

Munition 

HC 

WP 

Buildup  Time 

1.0 

0,5 

TABLE  14.  MUNITION  AVERAGE  BURN  TIME  (MIN) 


Gun 

105  Howitzer 

155  Howitzer 

Munition 

HC  WP 

HC  WP 

Average  Burn  Time 

3.0  0.0167 

4.0  0.0167 

TABLE  15.  RATE  OF 

FIRE  (ROUNDS/MIN)  VS  STABILITY  CATEGORY 

Gun 

105  Howitzer 

155  Howitzer 

Munition 

HC 

WP 

Stability  Catego*'y 

A B C D E F 

A B C D E 

i 
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CO 

Cl 

VO 

TO 

T1 

DO 

SO 

PO 

MO 

RO 

XO 

T2 

Y 

Z 

PSCTAB(7,9) 

T(5.4) 

C(5,4) 

B(5) 

G 5) 

H(5) 

R(5) 

CS(5,4,3) 


0(5,4) 


■ 

S(6.4) 


Y1 

-Y4 

R< 

2,2) 

G 

2,2 

F(2,2) 

H(2,2) 

A( 

2,2) 

B 

2,6) 

D 

2,6 

E( 

2,6 

E(2,6) 

W(2,2) 

SITE 

PSC(6) 


WLNGTH(5) 

P 

D 


Ceiling  - hundreds  of  feet 

Cloud  cover  - per  cent 

Visibility  - miles 

Temperature  - degrees  F 

Dew  Point  - degrees  F 

Wind  direction  - tens  of  degrees 

Wind  Speed  Knots 

Atmospheric  stability  category 

Mixing  depth  height  - meters 

Relative  humidity  - percent 

Total  distance  to  be  smoked  - meters 

Time  smoke  required  - minutes 

Average  roughness  element  - centimeters 

Roughness  length  - centimeters 

Table  of  stability  categories  (depending  upon  solar 
altitude  and  wind  speed) 

Table  of  trarsmittances  owing  to  water  vapor,  haze/fog, 
rain  and  smoke  for  0.55,1.06,  2.3,  3.8,  and  10.6  micrometers 
Table  of  smoke  concentration  values  for  fog  oil,  HC,  FS,  and 
WP  for  0.55,  1.06,  2.3,  3.8  and  10.6  micrometers 
Error  function  absorption  coefficients 
Scale  'leight  for  Mie  scattering  ' 

Haze  and  fog  attenuation  coefficients 
Rain  attenuation  coefficients 

Table  of  coefficients  used  to  calculate  smoke  concentrations 
using  the  calculated  transmittance  values  for  0.55,  1.06, 
2.3,  3.8,  and  10.6  micrometers 

Smoke  source  strength  values  for  fog  oil,  HC,  FS,  and  WP 

for  0.55,  1.06,  2.3,  3.8,  and  10.6  micrometers 

Coefficients  to  compute  ay  continuous  source 

Coefficients  of  the  roughness  correction  factor  used  in 

calculating  az  for  the  various  roughness  lengths 

Yield  factors  for  fog  oil,  HC,  FS,  and  WP 

Total  number  of  rounds  required  to  maintain  smoke  screen 

Number  of  guns 

Number  of  rounds  per  gun 

Unit  (per  gun)  source  strength 

Smoke  build-up  time 

Munition  average  burn  time 

Rate  cf  fire  vs  stability  category  for  105  Howitzer 

Rate  of  fire  vs  stability  category  for  155  Howitzer 

Total  time  for  munition  replenishment 

Met  observation  station  identifier 

Stability  category  indicator 

Wavelength  indicator 

Precipitation  indicator 

Demo  indicator 
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0 OaLi.oLmIE  iKAHbwIT  i ancE  OmIi'^o  TO  A T 1 LijUA  T 1 ON  BY  HaZl  aNU  FOG. 

1»-(P  .i_u.  Nu)  ou  To  3900 
|(j*2)=i.O 

oo  to  moo  ; 

3900  oOhI1i>(Ul 

lElv/l)  .uT.  o(J))  Go  TO  4uOu 
U0=H4 

l1=ho 

L2=U0 

1.3=0. b*lLi+L2) 

1.4=L2-L1 

i.b=i)  ,i;obb  ^bi*L4 

|3=ii.b»i.4*(»-Nt(L3+Lb)+FlML(L3-Lb)  ) I 

|4=tLXp(-HlJ)*T3)  ‘ 1 

L0=H3-rm 
ul=n4 

l2=M4+lo 
1.3=0. b*lLi+L2) 
u4=L2-i.1 
1.3=0.  «;oob7bl*L4 

|b=o.b*L4*(FNL)lL3+Lb)+FNL»(L3-Lb)  ) 
rb=tAP(-0.l2d*Tb) 

r(j»2)=|4*To 
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I 


OU  lu  4lOU 
40UU  uOu  I iUUt. 

LOihi 

1.1-MU 

i.i:=L0 

l.^=Utb^lLi+L«i) 

L4::L2-Li 

i.b-u«toob/bJ,*L*4 

1 7 = UtO^L4*(H-(U(Ub+Lb)+Kl>iu(LJ-Lb)  ) 

I t 7) 

4iUu  I i wOi. 

iH  tC-vji*  Yub)  bU  TO  HiiOO 
wKi  1 1.  (u»  bdUUO  ) r(j»^) 

b ChLCULaIL  li<A.4bivin  I AhCt  Owllxb  TO  ATTLuOAflON  bY  KAlN. 


4^Uu 

bON  1 l|\(Ut 

It"  ( M ,1.4  • Ybb  ) 

bo  TO 

44UU 

Abuo 

ObM 1 iHUb 

1 ( . 
oO  10  4auu 

*+4uw 

4ON  1 iijob 
it-(VU  .of*  4iO»  ) 

bO  To 

430U 

1 (jO)-i_At-'(-hb+K(o)  ) 

40UU 

4O14  ( i WOb 

it*  (iJ  .bui  Tub) 

bO  TO 

•P- 

c 

c 

c 

MKi lt(u»bVUuO) 

1 ( J • 0 ) 

C bMLbvjLAlt  I MAi'^bwIl  1 ANbt.  OWlUO  TO  ATTLiTUATION  bY  SMOkK, 
40UU  bO|4  I X|>iUt_ 

I tj»4)iu»U<d/(T  (»J»l)*T(J»£;)*f  (J»b)  ) 

•L.b»  !•)  00  To  47ou 

I U»4)-i., 

•+  7Uy  bUH  I INUL 

ll'CiJ  (bu*  Yc-b)  bo  |0  ‘♦oOu 
MbirL(u>4uUJU)  r(om) 
mKI  I L (o>  lU<iuO  ) 

b bAUbUbrtft  ji*iOkL  CbhbLHTKATION, 

40Uu  bUMlihOt 

.HL.  1.)  00  TO  bUUO 

l>0  UVu'j  1=1  » 4 
4 ( J» 1 ) -U • 

4':#Uu  bONllNOb 

bO  10  b^lUO 
bUuo  bO|\lll|.tUL 

[ b=ALOo ( T ( J»  4 ) ) 

iy=rb+io 
UO  bluu  K=l»4 

b ( J»K)=bb( J»K»l)+bb( J»K»^)*Ttt+Cb( J»K»3)*Ty 
blUU  bOi'^llNOL 
O^UO  bOi^TlNOL 

IKO  ,bO.  Yt-b)  bO  TO  b300 
wKl rL(o»4iU0U)  b(J»l) 


i 
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nvKi  ru(ui4«:uuO ) 

mKI |L(u»4jU00 ) L(J»^) 

WKl  IL(u>4<4U00  ) L(J»H) 
t^JUu  LUi>|I1|mUL 

3HUU  C*J|^T  li'jUt 

iH  (IJ  (CU*  fLb)  OU  TO  bjUO 
MKlILlotlUlUO) 

DbUo  cUi'l  I IimUc 

kl  r iJKi< 

C^^wKwaT  a I t-Mt.l'J  I b t 

lOuuU  f-Urti'-iA]  ( IHI ) 
lOioU  (-UKi'^Al  I IHU  ) 
iUiuU  rOKKiAl  UH  ) 
lOooO  hUHi'iA  I Ih  lU  . U ) 

bOUuU  I- UKi'iA  M <.t>M  A 1 1'lubKHLKlC  Ot-lACb  ANL)» 

^ j4n  bMUKL  tONCLNl kaT lOu  LALCULATlONb: ) 

biuoU  »-Ui<i'tA  I I jln  bUAi>j|  aAuOL  I (J  lAKbcI  - Mt-TLKb) 
b<iUuO  h UKi'iA  n jlri  aNOlC.  bh  bibur  I 1Aa(jc.I  “ UtG) 
jAUjU  KOKi'^AllAyn  bLANT  AAUOt.  TO  TAKOuT 

jOUuU  AUK'tAMAyn  ANoLL  OA  blohT  TU  TakGLT 

bUUuU  t-OKMAMnyn  pKt-O  IPl  I aOLL.  rtATLK 

JUboU  r- UK''"A  I (H^n  AMuUNT  Oh  wAltK  VAPOH  IN  PATl 

b/UuU  hUAi»iA|(H7n  I KANbMl  T f ANtt.  OWINo  10  ATTENI 


Ou  TO  bbUo 


bLANT  AANOt  TO  TAKouT  - MlTLHS  = *hht.^) 

ANoLL  OA  blohT  To  TakGLT  “ ULo  - 

pKt-ClPl  I AbLt.  rtATLK  - Cm/km  = »h8,^) 

AMOUNT  Oh  wAltK  VAPOK  IN  PATH  - CM  = » h B . ) 

IKANbMlTfANCt  OWINO  lU  ATTENUATION  oY t » 
l4fi»»A|LK  VAPuK  = rho.i;) 


bCOoU  hOKi*iA  I Ih7a»  i4hliA^L/A0G  = *(-0.2) 

byUoO  hOK-'IAMA7A»i4ht<Alix  - »ht)t2) 

hOOuU  h OkmA  1 I a7a  HHhlbMOKL  - »|-0.<d) 

4lUwO  hOKMAl  ( /A»4cHbi''ioKL  CONClN  1 K A T 1 ON  : (-UG  C 
H^UUD  hOKi*iA  I (i>VA»tUhML  - OM/bU  M = »Fo,2) 

40UOU  h Oki»IA  I I ^^A  » 20hh  b - oi^i/bO  M = »KO,2) 

‘♦‘♦UuU  hOKwAI  l^iyArcOhlwP  - OW/bu  M = »Kb,2) 

4oOoU  h OkmA I ( 4X » Ah » ion  mIcKOMo foKb; ) 
tNo 


FUG  OIL  ~ GM/bO  M = rhb.t) 


cjUuKOullNt. 
lwrtOt.K  PO 

V.OMi'ltiN  /^wl^/  u »CU » cl’ » '/U  * P » T U » T i » ijO  » bO » Y * PO » KU  * 

£1  jlbi4J*C(bi4)»AU»'iitb»4) 

/OUTPUT/  WLNoTH(b) 
ulrttiMbluN  Ato)»b(u»4) 

uAfA  m/u  •4»U.b*:»  144  0.076/ 

uAfA  b/U»llt  »U.lo»0.11ii»U.oyU,u.Uo09»0.06i8» 

<i  i.oo»o.yb»u«yJi»u.aoy»o.6yb»u.  763* 

3 o»3bt-04  »o«  b^E-04»  y . 06L"U4  > 1 ,3bE“03»  1 .ybt--u3»  1 » 361.-0 •>» 

4 0*Olb>0.7b»0.7l6»o.6OOr0. 604  r 0 * 672/ 


ijAfA  PI  /3.i4j.by£;ob4/ 
pAfA  KCb/oHYCb/ 

I i*lojPtlc.Kiu  UiFPUblON  Ai'^U  bMUKt.  buUKCt.  bTKtNbTH  CALCULATIONS. 
..Pi  ILluf  lOOOU) 
lH(iJ  .Lu»  Tub)  60  TO  bUUO 
i I L ( u»  bOOuU  ) 
wr\i  IHorlUilO'j; 
oOOU  oUl4  I il'jOL 

m61  I L.  ( 3»  blOUU  ) 
i\LmI  » ( u»  l*iOUU  ) AO 


.ft'  1 I L I J » blbou  ) 
ut-Ml)  I L>»  i<iUOU  ) MU 
«K1 IL(b»b20u0) 
KtAiilbf  l^iOOO  ) £.0 
k»K1  lt(3»3<ibU0) 
i\Lmm  ( b»  iiiuOO  ) AU 

.Lu*  Tub)  60 
Mt<l  I L (u#  bjUuG  ) AO 
mKx I L ( U . bObuO ) MU 
««Ki  lL(u>b4UuO  ) ^.U 
mM  I L ( o»  b4bU0  ) Au 


lO  bxuu 


blUo  oOi4 1 iu06 

C oiKFOblON  OAUOULATlUNb  FOK  COlMTlNUOUb  bOUhCL. 
Ai--1 .tH+l.iy^ALOolut Y ) 


i.=lij  . ♦♦aI 
ul=ALo6 ( ^ ) 


o«i=LSl*bl 


u3-o2*ol 

D4=M3*di 

ab=M44ui 

ub=u.^4^Oobo6y+U.^y404y£^6b♦bl-0. 23721 3914*82 
l>7=u»  Ibb34yb04 *63-0. 03201  b723*d‘++2.  lblbbt.-03*bb 
Ul^ttb-fd/ 

0l=tXp(Ul) 

oO=-l.^yb£i8byuy-1.0u6lOb7b4*dl  + 1.4db094b8b*b2 
u 7--U . 774l3o72b*63+U • lb6bb93bb*b4-0 .010823351 *bb 
U2=6b-fu  7 

u2=tAP(u2)-0.22b 


\ 

I 


I 


i 


11- .vjF»  y.yvyyyy)  Gu  to  b^iuo 
ijb=b,77^o7t.-04+«;.51y‘lJ£-0b*dl  + b.7lUHiE“Ub*G2 
u /^-b.HoOU^L-OO^LiO  + i # 042 1L“0 7*04  + 2.  bbl 31  t-0tt*bb 
oU  10  o.:<:u 
o2Uu  OOtNillUOL 

IH  (2  ,ol  . 40.  ) GO  To  b<:lu 

Dt>=-li,oOi34yoi  + 2.14G242bl4*bl-U.  ibb210bl7*U2 
L>7=7.ooo0^b.-O3*b3-I.473b3c-U4*d4+1.  lo^bbe-Ob*Bb 
oO  10  b£;2u 
O^lo  OOu I iuOL 

uD  = lluo.3ubbtio-iU0.  j4ytJo3o*bl  + 2.4244y42b6*LJ2 
b /=-0  . ui4bdH77b*bO+4 . 34bl  7t.-ob*b4-4  • o9bbbE-0ti*Hb 
oOH  I iuut- 
bO-Mb+L)  / 

bb-U. 0007/00  ljy  + i.oy^bl4  /bo  *01-1 . b730obB3b*G2 
b / = U . /bO  7y*b.)-0 . 140b<;oyu4*b4  + y .blb21L-U3*Mb 

o4-bb+b / 
o4it.AH(u4)-i 

!>-(/:  .bl.  10.)  OO  10  b2bU 
bl=rtLoO lUi*A0**o2*l/(i+U0*XU**04) ) 

OO  10  o<;4u 
0<^ClO  OOn  I ll.UL. 

bi-ALoo  <01*aU*»U2*  ( 1 + 1/  (l»J*XO*  *04)  ) ) 

J^‘*U  OON  r ilNOL 

b<i-olKO»l)*A0**b(HU»<i)/(i  + blP0»b)*AU**b<PU»4)) 
jc=nl+b2 

oi-A (ho ) ♦ao**o  .y 

l)'(l>  .(.o.  Tub)  OO  |0  0200 
mHI (E (u. 1o2oU ) 
mKI lE(u.boOOO) 
wKl I t (U.1020U) 
wKl I t (b.bbOoO)  bl 
mKX lE(o.bob00 ) b2 
0^00  OOwllUOE 

0 bi'',OM.  bOUKOE  bFKElMolH  CALCULA  1 l0|N(b  H OK  EUNTINUOUb  SoUkCE. 
H^i^AboiMO-OU  )»  (Hl/lbO.  ) 
ll-(yU  .kE.  0.)  oo  To  b3oo 
bO-1  .u 

odou  ^OhIIwOE 

j3=so*o .bib 

«U=y2*00*b0KT  (h1  ) /boKT  ( 2 ) ♦EXH  (-0  . b*  ( (20“H0  ) /S2  > **2  ) 

Ki:=>>uKr  (b.3124/  (b.7o*Si|>i(M2)*blN(  A2)+0.49*COb(A2)*COb(  A2)  ) ) 
r«:=U.yjj7+(u.U3by*Ko)-(7b-04*hU*K0)  + (b.llE-0b*H0*K0*K0) 
Tb=1.07/b+(0.uyob6*K0)-(l.bE-0-i*K0*K0)  + (1.56E-0b*M0*»<0*H0) 
wO=uO*K2 

wi=oo/r2 

o2=oO/Y3 
uO  0400  l=l»b 
o(l>l)=o(l*l) *OU 
Vi<(l»2)=o(i»2)*ol 


1 


i 

j 
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5 


if-dJ  .uu.  YLb)  t.0  TO  b<+UU 
mKI  IL(u»1Uc:UU) 
mKx  [LiorbS'OUO)  W|_ImGTH(1) 
wHilLlu>lu<^UO) 

MKi It (u»boUuO)  U(i»i) 
mKi I L (u* aubUO ) U(i*^) 

•vKx  It  ( u»  b/UUO  ) vj(irb) 

OHUU  I ilMUU 

C uiFl-UbxuiM  L/»LtOLAT  lUiJb  l-OH  UOAbx-lNbT  AimTmNLOUS  bOUKCF, 
bi-^./j. ♦bl 

* / -ii » ♦:5<i 

ll"  (u  ,Lu»  iTtb)  Oo  10  bbUu 
mKI  IL(utlU<iOO) 

AK|lt(u»oiUUO) 

X 1 1.  ( u I lO^uU  ) 

X I c.  ( u » bbUuU  ) bl 
.•AX  lL(U»bbbOU)  b<: 

:j boo  cOki  I XuUL 

c bi'iOM.  bOUrxcu  ;jTAtUoIn  CaLCOLaT  l0i<b  H UK  UUAbl-INbTANTAWtOOS  bOUKUt 

v«o-tbl*3<i^Hx)/LAK(-lZ0*^iJ-M0»hU)/(i;*:j2*b2)) 

I 4^  4u'i-*-  ( 0 .U  / /‘♦♦Ko  ) - ( i .ot.-0b*l<U*K0  ) + ( 1 • 73t-0b»Ku  ♦KU*K0  ) 

jO-l^*U/Y4J*XU 
t>ooo  i-l»o 
ui  ( i f 4 ) -c  ( X » 4 ) 

xHli)  .L.U*  Yi_b)  OU  |0  boUU 

mKX lt(u»lU20U) 

mKx  It.Tutb'^UUU  ) Mti^tV>TH(l) 

mKX  It.(u#llj200l 

•VK  X 1 1 ( u * bosuU  ) u(l»4) 

3UUO  V.UI4  I iUUL. 

it"  In  .tv<*  Yt.b)  OU  lO  b7uo 


••KX  I L (U»  luiuo  ) 
o7Uo  cUwliuUL 
rtt  r UKN 

t*h^Ki<iHT  b I A I b « 

iOuoO  l-Oiu'iA  I ( XtU  ) 
lOloU  HUrti''iA|  UNO) 

iOcOU  rUKi'-iAl  (XH  ) 

12UuO  I-Uki''IAI  d-iU.U  ) 

bOuuU  eUKMAI l4A»boHA|MUbPMEKlC  ulFKUblOu  AND  bMOKE  SOUKCE» 
i>  blKLNuTh  CALLULATiOWb:  ) 

blUuO  l-UKiviAno2ri  lOTAL  UlblANCL  TO  bt  bMOKLD  - M) 
blboO  FUrtMATIblM  KtLtAbt  hEIOhT  OF  bMoKt  bOUKCE) 
o20u0  FOKMAKblH  MEAN  HtloHJ  OF  TAKOET  - MLTEKS) 
b^DuO  F Oki'^A  I ( bid  UlKtCTiON  OF  LlNt  oF  blGMl—OEG) 

bOUUO  FOkwATC  7X»27mTo1  At  UlbTANtt  TO  bt  bMOKEL>»  iOX » ilH-  METERS  = »FH.2) 
obooU  t-OKWAl  (7X»3oMKLLLAbc  HEIGHT  OF  bMOKE  SOURCE  (AGL)» 

*;  i2H  - meters  = ff-o,id) 

oAOoO  hUKMAI (7X»2iHMtAN  HtiGHT  OF  TARGE T » IGX » 1 IH-  METERS  = »FB.2) 


I 


i>UUKUu  I NrtlKH 
il'IltOLK  Hu 

cUiV|iviUu  /KiviK/  u » CU  H,  1 » '/U  » H » I U r Ti  » UO  » bO  » Y »P0  » KO  » 

I ( L) » *+ ) » C ( i5  » 4 > » AO  » ^ J tJ  » 4 ) 

cUiviMUh  /ournul/  wlNoTH(O) 

L) i McNs  1 ui'^ 

jllxitNiluN  lOi.  (t»0)  »L(2»0)  »|-  fhiliiP 

JA|  A ( •I  = l»*i)/it3.y»i7i7.J5»4b,tt»7oyb,2/ 

u^fA  I iM(lfJ)  * 1-1 1 c.)  / i t »u»a#l»  »U*b/ 

uAfA  ( fJ-1  fid)  »I=l»<;)/b,  »u«u1d7»4.  »U.01o7/ 

ijA(a  l(oi(i»J)«o::ifti)*I  = i»<^)7b>»4,»b>»2>*li»li>o*»0«*t}*>4,>ltb*l« 
u^lA  ( Il(  1»J)  f 1-1 1 ^ ) / ^ , *«;,  »0»!3»U»333» 

i U • r u . > 3 . f 2 , > 0 • b f 0 t 333/ 

bA  I A Yt-b/jHYbb/ 

I lOiJb  uAPbl'Ib  1 I UKub  * 
aKI I t (u» lUUuO ) 
mbx lUlurbJUjO) 

.•b  i ( t ( u > Iu2u0  ) 
aH  X I (.  ( j « bUbuU  ) 
i\t.Ai  <(Li>x<^UOij)  !<:: 

L^O  oUOU  A — 1 » b 

iJ'tb  .I-V..  1)  bu  lu  uOOU 
xt-t.).  Lw«  YLb)  bu  TO  byuo 
oOuu  bUu  I ii^uL. 

.t-U.  i.)  uO  TU  o7b0 
b bMLbbl-HiL  uuMotb  Oh  bUNi  KCuDlKtU, 

b(ir2)-u(b»b)/tl(l*i::) 

bl4i»2)-v.(K»b)/h(c;»2) 
ub  n2uu  l-i*2 
uO  blou 

bU-Ai|\i|  (b(  1 » J)  ) 
bl-u ( X » b ) “bU 

il"  (bl  ttb»  U*)  bU  Tb  bluU 
b ( 1 • 0 ) ibO+i 
olUo  bbijIii^Uu 
O^UO  bbN  I il'JUL 

b bMUbbUAlt  IuTaL  Hail  hUl<  KtPLLWlbHMLNT . 
bU  o4uU  1-1 
uO  o30U  J=lr2 
w ( X 'wl)  = I 2+A  ( I » J ) -t)  ( i » J) 
u3UU  LUi'xTInUL 
b40u  bUl^4li|^JUL 
b bMLbULAlL  KaTL  UF  FlKL. 

UU  obuu 

F HtU)  = (*<af  JUUKJiPo)  )4i 

Fl2»J)=(W(2»J)»L(J*H0))+l 
obUu  bbiXlihUL 
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C CmLCULaIU  IoTmL  NOMbtK  oF  KOUNUb  KtUUlKEU, 
ijU  n7uU  1 — 

JU  nbUU 

Ki=Aiix,r(K(I»J)  ) 

Ki:-K(  i f j)-Kl 

.c.Ui  U*)  bU  TU  bbuO 
K< 1 » J)-k1+1 
bbUu  CUiMT 

b7Uu  C'Ji'l  I iNUu 

(l)  (bu*  Yt-b)  Ou  lO  b7Jo 
aH  1 I k.  ( D > bHbuO  ) wLn07H(K) 
mHI  It.(u»10^u0) 
o7oyj  bbi<4liNUl. 

mKI  lt.(o»biUU0) 

|L(b>blbU0)  XU 

wKl  lt.(D»O^U00) 

mKI  I L (b» lU^uO ) 
wKi lb(o»b^bU0) 

wKi  li.(u»obUuO)  b(l»l)»LJi(l»KO)»K(i»l) 
wKl  lL(u>bobuO)  b(i;»i)  »t.  (i»PU)  »K(2»1) 
nKilL(u»lU^UO) 
mKI IL(o»b40uU) 

wKi lt(u»bibU0)  XQ 
mKx  I L io»  OcilUuO  ) 
wK i I L ( u » 1 U^UO  ) 
wKi IL(o»b^buO ) 

wKllt(orOOOuO)  b(l»^)  ibH^tHu)  »K(1  fid) 

«H1  It,(o»bbbu0)  b(^»^)  »E(,d»Ku)  »K(2»ii) 
bO  lu  o7yo 
b7bu  bUwl iNOb 

iH  (L)  tbu*  YbS)  bu  |0  b7dU 
wKi  lt.(o#b4bU0  ) wLnG7H(K) 
mKI lt(b»lU^UU) 
b7ou  cONTlNUt 

mKI TLlbfbbUUO) 
b^Sb  UUNriNOL 

mKI  rt.(brlO<iUO) 

btiUU  uUiMlii>jUL 

uyuu  cbiXfiNUb 

mKI IL(bf lUUUO) 

KtrUKN 

C*»-wKmaT  blArtMLNTS. 
lOUuO  FUKMAI (iHi) 
lOiuU  (-OKMAftiHO) 
i02O0  hOKMATdH  ) 

120u0  FOkmAT (FiOtO) 

bOUuO  FOKMAT(4Xfi2bHMUNlTlUNS  LKPEiNUlTURtS: ) 
bObUO  F'OKMATlblH  TlMt  bWOKE  KtUOlKtU  - MINUTES?) 
biUUO  FOKmAI (4X» JlHVOLUME  OF  FiKt  - HC  SMOKESCREEN) 


AQ 


■ 


oiiJUO  KUKi'^iaH /X»iioMbLKLt.N  LEN(»TH  - MtTtKb  = »F8,2) 

o2UuU  KUkmAT (7X»26HbCKLLN  UUKATION  - MINUTLS  = rF8.2) 
o2t)u0  i-OKMAl  lit>X»HHOUNb»t>X»lUHKOUNUb/MIW»bXil2HT0TAL  KOUNUb) 

uOuoU  t-UKi*lAT  t/X»^MlObOA»Ffc>«2»luX»Fb»^»llX»F6t2) 
uODuU  t-OKi'irtl  ( 7X » iMlbtjr  ^X»F^>»2»  luX  »Fb»2»  IIX  »F6»2) 
u‘*uuO  hOKi'iAl  li*Xt^lHVULUNlL  OF  F IKt  - »(H  SMOKtSCRtEN ) 
oAooU  FOki“iAT  (4*X*Ah»  JLoM  MlcROMtTLKb;  ) 

obUUO  FO«MMl(‘4X»4VHbi''iUKL  nOT  KLuOlKtO  DUE  TO  ATMOSPHERIC  CONDITlONb.) 
tNu 
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BASIC 

GLOSSARY  OF  MNEMONICS 


CO 

Cl 

VO 

TO 

T1 

DO 

SO 

PO 

MO 

RO 

XO 

T2 

Y 

Z 

P(7.9) 

T(5.4) 

C(5,4) 

3(5) 

G(5) 

H(5) 

R(5) 

0(4,3) 


Q(5,4) 

A(6) 

S(6.4) 

Y1-Y4 

R(2,2) 

G(2,2) 

F(2,2) 

H(2.2) 

A(2.2) 

B(2.6) 

0(2.6) 

E(2.6) 

W(2.2) 

I$(3) 

0$(6) 

A$(21) 

P 

0 


Ceiling  - hundreds  of  feet 

Cloud  cover  - per  cent 

Visibility  - miles 

Temperature  - degrees  F 

Oew  Point  - degrees  F 

Wind  direction  - tens  of  degrees 

Wind  Speed  - knots 

Atmospheric  stability  category 

Mixing  depth  height  - meters 

Relative  humidity  - percent 

Total  distance  to  be  smoked  - meters 

Time  smoke  required  - minutes 

Average  roughness  element  - centimeters 

Roughness  length  - centimeters 

Table  of  stability  categories  (depending  upon  solar 
altitude  and  wind  speed) 

Table  of  transmi ttances  owing  to  water  vapor,  haze/for, 

rair.  and  smoke  for  0.55,  1.06,  2.3,  3.8,  and  10.6  micrometers 

Table  of  smoke  concentration  values  for  fog  oil,  HC,  FS,  and 

WP  for  0.55,  1.06,  2.3,  3.8  and  10.6  micrometers 

Error  function  absorption  coefficients 

Scale  height  for  Mie  scattering 

Haze  and  fog  attenuation  coefficients 

Rain  attenuation  coefficients 

Table  of  coefficients  used  to  calculate  smoke  concentrations 
using  the  calculated  transmittance  values  of  0.55,  1.06 
2.3,  3.8,  and  10.6  micrometers 

Smoke  source  strength  values  for  fog  oil,  HC,  FS,  and  WP 

for  0.55,  1.06,  2.3,  3.8  and  10.6  micrometers 

Coefficients  to  compute  oy  continuous  source 

Coefficients  of  the  roughness  correction  factor  used  in 

calculating  az  for  the  various  roughness  lengths 

Yield  factors  for  fog  oil,  HC,  FS,  and  WP 

Total  number  of  rounds  required  to  maintain  smoke  screen 

Number  of  guns 

Number  for  rounds  per  gun 

Unit  (per  gun)  source  strength 

Smoke  build-up  time 

Munition  average  burn  time 

Rate  of  fire  vs  stability  category  for  105  Howitzer 

Rate  of  fire  vs  stability  category  for  155  Howitzer 

Total  time  for  munition  reple.iishment 

Met  observation  station  identifier 

Stability  category  indicator 

Wavelength  indicator 

Percipi tation  indicator 

Demo  indicator 


10  con  D,G0,C1,V0,P,T0,T1,D0,S0,Y,P0,R0,T[5,4],C[5,4],X0,Q[5,4] 

20  mOM  lOVIK:  IffiTEOROLOGICAL  INPUTS  AND  METEOROLOGICAL  CALCULATIONS. 
30  DIM  r.:f[5],P[7,9],Q$[6] 

40  FIXED  2 
50  PRINT 
60  PRINT 

70  PRINT  "KV/IK  SMOKE  PROGRAM" 

80  PRINT 
90  PRINT 

100  DISP  "IS  THIS  A DEMO  - 1=YES  0=N0" ; 

110  INPUT  D 

120  DISP  "MET  SITE  ID"; 

130  INPUT  I;; 

140  DISP  "LATITUDE  OF  MET  SITE  - DEG"; 

150  INPUT  T.0 

160  DISP  "LONGITUDE  OF  MET  SITE  - DEG"; 


170 

TNPU'n 

LI 

180 

DISP  " 

'ALTITUDE  OF  MET  SITE-KILOMETERS"; 

190 

INPUT 

70 

200 

DISP  " 

'cNJLIAN  DATE  OF  MET 

OBSERVATION"; 

210 

INPUT 

JO 

220 

DISP  " 

'ZULU 

TIME  O]'’  MET  OBSERVATION-HR"; 

230 

INPUT 

HO 

240 

PRINT 

It 

MET  SITE:" 

250 

PRINT 

260 

PRINT 

tt 

ID 

= ":I 

270 

PRINT 

(1 

IJVTITUDE 

- DEG 

= ";L0 

230 

PRINT 

tt 

LONGITUDE 

- DEG 

= ";L1 

290 

PRINT 

It 

ALTITUDE 

- KM 

= ";Z0 

300 

PRINT 

310 

PRINT 

tt 

JULIAN  DATE 

; - DAY 

= ";J0 

320 

PRINT 

It 

ZULU  TIME 

- HOUR 

= ";H0 

330 

PRINT 

340 

PRINT 

350 

DISP  " 

'CEILING  - HUNTiREDS 

OP  FEET' 

It  • 

> 

360  INJUT  CO 

370  C0=(C0*100)*0.3048 

380  DISP  "CLOUD  COVER  - PERCENT"; 

390  INPUT  C1 

400  DISP  "VISIBILITY  - MILES"; 

410  INPUT  VO 
420  VO=VO*1.61 

430  DISP  "PRECIPITATION  - 1 =YES  0=N0" ; 
440  INPUT  P 

450  DISP  "TEMPERATURE  - DEG  P"; 

460  INPUT  TO 

470  TO=(5/9)*(TO-32) 

480  DISP  "DEW  POINT  - DEG  F" ; 

490  INPUT  T1 

500  T1=(5/9)*(T1-32) 
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510  DISP  "WIND  DIrtEOTION  - TJINS  OF  DEOS"; 
520  INPUT  DO 
530  DO=D0*1O 

540  DISP  "\,'IND  SPEED  - KNOTS"; 

550  INPUT  SO 

560  DISP  "AVE  ROUGHNESS  ELEMENT  - CM"; 

570  INPUT  Y 

580  PRINT  " I4ETE0R0L0GICAL  INPUTS:" 
590  PRINT 


600 

PRINT  " 

CEILING 

- 

METERS 

= 

II 

CO 

610 

PRINT  " 

CLOUD  COVER 

- 

PERCENT 

= 

II 

01 

620 

PRINT  " 

VISIBILITY 

- 

KILOMETERS 

s 

II 

VO 

630 

PRINT  " 

PREOIPITATION 

= 

II 

p 

640 

PRINT  " 

TEMPERATURE 

- 

DEG  C 

ss 

It 

TO 

650 

PRINT  " 

DEWPOINT 

- 

DEG  C 

= 

11 

T1 

660 

PRINT  " 

WIND  DIIWCTION 

- 

DEG 

= 

It 

DO 

670 

PRIKroi  ,t 

WIND  SPEED 

- 

KNOTS 

= 

rr 

SO 

680 

PRINT  " 

AVE  ROUGHNESS  ELEJP-JNT 

- 

CM 

= 

II 

Y 

69C  PRINT 
700  PRINT 
710  FOR  J=1  TO  9 
720  FOR  1=1  TO  7 
730  FJ-LAD  P[T,J] 

740  Ni^YT  I 
750  NT-ECT  J 
760  ra^AD  ''''Til 
770  IF  01 #100  TUFF  820 
7B0  IF  C0>2133.6042  THEN  820 
700  11=0 
000  12=0 
810  GOTO  1540 

820  7KM  CALCUI.ATli  ANGULAi^  FRACTION  OF  A YTIAR  FOR  A GIVEN  JULIAN  DATS  (AO). 
850  R9=PI/180 
840  D9= 180/FI 
850  L0=T,0*R9 

860  AO=((JO-1 )*360)/365.242 

870  rffiM  CALCULATE  SOLAR  DECLINATION  ANGM  (A4). 

880  A1=A0*U9 
890  A2=279.934B+A0 

900  A2=A2+(1.914827*GIN(A1 ))-(0.079525*G0S(A1 )) 

910  A2=.A2+(0.019938*SIM(2*A1  ))-(0.00162*COS(2*AO) 

920  A2«A2*H9 

950  A5«23.4433*R9 

940  A4=SIN(A3)*SIN(A2) 

950  A4=ATN(A4/SQR( 1 - AA*A4+1E-99) ) 

I 960  REM  CALCULATE  THE  TIME  OP  MERIDIAN  fasSAGE  - TRUE  SOLAR  NOON  (A5). 

I 970  A5=12+(0.12357*SIN(A1 ))-(0.004289*C0S(A1 )) 

980  A5=A5+(0.153809*SIN(2*A1 ))+(0.060783*C0S(2*A1 )) 
i 990  REM  CALCULATE  SOLAR  HOUR  ANGLE  (A6). 

I 1000  A6=15*(H0-A5)-L1 
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1010  A6=A6*H9 

1020  ret:  calculate  solar  altitude  (A7). 

1030  A7=SIN(L0)*STN(A4)+C0S(L0)*C0S(A4)*C0S(A6) 

1 040  A7=ATN( A7/S0R( 1 -A7*A7+1 E-99 ) ) 

1050  REF  CALCUMTE  TIME  OE  SUNRISE  AND  ‘’UNSET  (BO,Bl). 
1060  AR=-1 .76459*Z0*0. 40705 
1070  A8=Aa*R0 

1080  A9=(SIN(A8)-(SIN(L0)*SIN(A4)'))/(G0S(L0)*C0S(A4)) 

1 090  A9=ATN  ( S':R ( 1 -A9*  AO ) / ( A9+1 E-09 ) )+2* ATN1 E+90*  ( A 9<0 ) 

1100  A9=A9*D9 

1110  A9=A0*( 24/360) 

1120  !^0=fLl/lS)-‘-A5_Aq 
11^0  B1  = (L1/15)+''S+A0 
1140  IE  81  <-  24  TIUDN  IK^O 
1160  131=81 -'5/ 

1160  A7=A'7*D9 

1170  REM  OAT-OULATE  INSOLATION  CLASS  TTUKPER. 

1180  12=0 

1100  IE  A7  <=  60  THEN  1220 
1200  12=4 
1210  GOTO  1300 
1220  IE  A7  <=  33  TUE'N  1250 
1230  12=3 
1240  GOTO  1300 
1250  IE  A7  <=  15  then  1280 
1260  12=2 
1270  GOTO  1^00 
1280  IE  A7  <=  0 TH'O'I  1400 
1200  12=1 

1300  TWi  CALCUT/T]’:  NST  ILADTATIOM  liOEX  EOR  DAYTTTU!;. 
1310  13=0 

1320  IE  01 >50  THEN  1350 
1330  13=12 
1340  GOTO  1430 

1330  IE  CO  >=  2133.6042  THEN  1380 

1360  13=12-2 

1370  GOTO  1430 

1380  IP  CO  >=  4876.8096  THEN  1410 
1300  13=12-1 
1400  GOTO  1430 
1410  IE  Cl, f 100  THEN  1430 
1420  13=12-1 
1430  IP  13/0  'THEN  U50 
1440  13=12 

1450  IP  I3>1  THEN  1470 
1460  13=1 
1470  11=13 
1480  GOA'o  1540 

1490  'lETI  CALCULATl'.  HET  RADIATION  INDEX  FOR  NIGHTTIME. 
1300  IF  Cl >40  'CHEN  1530 
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1510  I1=-2 
1520  CrOTO  1540 
1530  I1=-1 

1540  REH  CALnUMTE  PASC'UTLL  STABILITY  CATEGORY. 

1550  14=0 
1560  15=0 

1570  IP  11/^4  THEN  1590 
1580  14=1 

1590  IP  I1/;-‘3  THEN  1610 
1600  14=2 

1610  IP  I1t‘2  then  1630 
1620  14=3 

1630  TP  11 '^1  THEN  1650 
1640  14=4 

1630  TP  11 ^0  THEN  1670 


1660 

14=5 

1670 

IP  11 

,^-1  TTfSN  1690 

1 680 

14=6 

1690 

IP  11: 

/.-■-2  THEN  1710 

1700 

14=7 

1710 

IP  SO 

>=  2 

TJP'JN 

1730 

1720 

15=1 

1730 

IP  SO 

>=  4 

THEN 

1760 

1740 

15=2 

1750 

G0''’0 

1G50 

1760 

IP  SO 

>=  6 

TinCN 

1790 

1770 

15=3 

1780 

GOTO 

1950 

1790 

IP  30 

>=  7 

THEN 

1820 

1800 

15=4 

1810 

GOTO 

1950 

1820 

IP  SO 

>=  8 

'I'TDDN 

1850 

1830 

15=5 

1840 

GOTO 

1950 

1850 

IP  SO 

>=  16 

1 THEN 

■ 1880 

1860 

15=6 

1870 

GOTO 

1950 

1880 

IP  SO 

>=  11 

THEN 

1910 

1890 

15=7 

1900 

GOTO 

1950 

1910 

IP  SO 

>=  12 

! THEN 

1940 

1920 

15*8 

1930  GOTO  1950 

1940  15=9 
1950  P0=P[I4,I5] 

1960  REM  CALCULATE  MIXING  DEPTH  HEIGHT. 

1 970  M0= ( 6-PO )* 1 21 * ( T0-T1 )/6+ ( P0*0. 087* ( 30+0 . 5 ) )/ (1 2*8 . 237B-05* 3 . 809 ) 
1900  REM  CALCULATE  RELATIVE  HUMIDITY. 

1990  IP  TOO  THEN  2030 
2000  AO=9.5 
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1- 


( 

t 

2010  BO=265.5 

2020  GOTO  2050 

2030  AO=7.5 

2040  B0=237.3 

2050  XT'  T1>0  2090 

2060  A1=9.5 

2070  B1 =265.5 

! 2000  GOTO  2110 

! 2090  A1=7.5 

( ■ 2100  B1=237.3 

2110  EO=6.11*10*((AO*TO)/(BO+TO)) 

2120  E1=6.11*10*((A1*T1 )/(B1+T1  )) 

2130  R0=(El/i:0)*100 

2140  PRINT  " IPETROROLOGIOAL  CALCULATIONS:" 

2150  PRINT 

^160  PRINT  " PASOUTLI.  STAT3ILITY  CATEGORY  = " ;(i:i?[PO,POl 

'’170  PRINT  " RTiL-ATIVE  HUMILITY  = ";R0 

2180  PRINT 
2190  PRINT 

2200  LISP  "])ONE  - LINK  1" 

2210  REM  PAS  'UILL  STABILITY  CATEGORY  LATA. 

2220  LATA  1 , 1 , 2, 3, 4, 6, 6 
2230  LATA  1 , 2 , 2, 3, 4, 6, 6 
2240  LATA  1 , 2 , 3, 4, 4, 5, 6 
2250  LATA  2, , 3, 4, 4 , 5, 6 
2260  LATA  2,2, 3, 4,4, 4, 5 
2270  DATA  2, 3, 3, 4, 4, 4, 5 
2200  GATA  3, 3, 4, 4, 4, 4, 5 
2290  LATA  3, 3 , 4 , 4, 4, 4, 4 
2300  LATA  3, 4 , 4, 4 , 4, 4, 4 
2310  LATA  "ABCLNF" 

2520  HNL 


10  COM  l),OO,G1,VO,P,TO,T1,DO,SO,Y,PO,RO,T[5,4],C[5,4],XO,0[5,4] 

20  RJ'il'i  la/IK:  ATMOSPHERIC  OPTICS  MJD  SMOKE  CONCENTRATION  CALCULATIONS 
30  DIM  P[5l,cr5],n[5],R[5l,A::5[2l],D[4,3] 

40  FIX]<]D  2 

50  IP  D=1  THEN  80 

60  PRINT  " ATMOSPffiRIC  OPTICS  AND  SH0I3  CONCENTRATION  CALCULATIONS:'' 
70  PRINT 
80  FOR  1=1  TO  5 
90  READ  B[l],0[l] 

100  IRJXT  I 
110  V1=LOO(VO) 

120  V2=71*V1 
130  V3=V?*V1 

140  H[1 1 = 1 .55‘51-(0.9811*V1  )-(0.0197*V2)+(0.0041*V3) 

150  Hril=SX?(H[ll) 

160  Tir?'  =1 .5551-(o.9811*V1  )- (O. 01  97*V2 )+ (0.0041  *V3 ) 

170  Hr2l=EXP(Hr2ll 

180  Il[ 3'  =1 . 449i- ( 1 . 0044*V1  )- (0. 01 2*V2 )+  (0.0032* V5 ) 

190  Hr3'=EXP(ll[3l) 

200  Tir4'  =1.2304-(1 .0436*71  )+ ^0. 0099* Vf )_ (o. 001 6*V3 ) 

210  Hr4l=EXP'(N[4l) 

220  ir5'  = 1 .5176-(1 .7147*71 )+(0. 0001*72 )+(0. 0423*73) 

230  Hr5l=EXr(l![5l ) 

240  ur 1 ■ =1 . 3506- (0. 8825*71  )- (0. 0753*72 )+ (0.01 29*73 ) 

250  Rri'=i;XP(Rril) 

260  R[2l=1 .4098- (0.9865*71 )-(0. 014*72)+(2.3E-03*73) 

270  Ri  2l='V:''(R[2l) 

280  Rr3]  = 1 . 5497-^0.8696*71  )- (0. 1 004*72  U (0.0231 *73 ) 

290  Rr3l=EX"’(.R[3l ) 

300  'r4]  = 1 .5556- (0.0013*71  )-(0. 0775*72 )+(0. 0173*73) 

510  Rr4J=E:G’0^r4l  ) 

320  Rr5l=1.5923-(0. 9396*71 )-(0. 0627*7? )+(0. 0168*73 ) 

330  r[5]=EXP(  {[3]  ) 

340  TIO=0 

350  DISP  "SLANT  R.ANOE  TO  TARGET  - METERS"; 

360  INPUT  H3 

370  DISP  "ANGLE  OF  SIGHT  TO  TARGET  - DEG"; 

380  INPUT  S 

390  IP  D=1  THEN  430 

400  PRINT  " SLANT  RANGE  TO  TARGET  - METERS  = ";n3 

410  PRINT  " ANGLE  OP  SIGHT  TO  TARGET  - DEG  = ";S 

420  PRINT 

430  H3=H3/1000 

440  IP  S >=  0 THEN  460 

450  S=-S 

460  S-S*(PI/180) 

470  H1=SIN(S)*H3 
480  H?=C0S(S)*H3 
490  S-SIN(S) 

500  H4=0 


w 
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510  IP  0=0  THEN  540 
I 520  H4=1/S 

I 530  REM  CALCULATE  PRECIPITABLE  V/ATER. 

540  V/=0.4477+(0.0328*T1  )+(1.2E-03*T1*T1  )+(1.84E-05*T1*T1*T1) 

: 550  IP  D=1  THEN  570 

560  PRINT  " PRECIPITABLE  V/ATER  - CM/KM  = ";V/ 

i 570  REM  CALCULATE  .\M0UNT  OP  V/ATER  VAPOR  IN  PATH. 

580  UEP  PNA(A)=-EXP(-S*A/2) 

■Vqo  1,0=113 
600  L1=U0 
610  L2=L0 

620  L3=0.5*(L1+L2) 

630  L4=L?-L1 
640  1,5=0.2886751*14 
650  V/0=0.  5*L4*(PNA(L3+L5)+PNA(L5-L5)  ) 

660  V/1=V/*>/0 
670  IP  ri=1  T/IEN  690 

680  PRINT  " AI'IOUNT  OP  WATER  VAPOR  IN  PATH  - CM  = '';\/1 

690  REM  CALCULATE  TRANSMITTANCES  FOR  0.55,1.06,2.3,3.8,10.6  MICRON  WAVELENOTHS. 
700  FOR  ,1=1  TO  5 
710  REAL  AO 
720  IF  D=1  THEN  760 
730  PRINT 
740  PRINT  AO 
750  PRINT 

760  REM  (CALCULATE  TRANSMITTANCE  OV/ING  TO  ABSORPTION  BY  V/ATER  VAPOR. 

770  IF  J.f5  THEN  800 
780  T[J,1]=EXP(-0.0681*V/) 

790  GOTO  900 

000  DEE  FNB(3)=EXP(-R'‘2) 

81 0 L0=  ( Bf  J J -it.SQR(V/1  *PI ) /2 ) 

820  L1=H0 
830  L2=L0 

840  L3=0.5*(L1+L2) 

850  L4=L2-L1 
860  L5=0. 288675 1*L4 
870  T2=0. 5*L4*(PNB(L3+L5)+PNB(L3-L5) ) 

880  T[j,ll=(2/SQR(PI))*T2 
390  T[J,1 J=1-T[J,11 
900  IP  D=1  THEN  920 

910  PRINT  " TRANSMITTANCE  OV/ING  TO  ATTENUATION  BY:  V/ATER  VAPOR=  ":T[J,1] 

920  RET!  CALCULATE  TRANSMITTANCE  OV/ING  TO  ATTEVlUATION  BY  HAZE  AND  FOG. 

930  IP  P-0  THSN  960 
940  T[J,2]=1 
950  CrOTO  1250 

960  IP  VO  >=  G[J]  THEN  1170 
970  DEF  PNC(C)=EXP(+C*S*L0G(0.1/H[J] )) 

980  L0=H4 
990  L1-H0 
ioOO  L2-L0 
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1010  ■'.3=0.5*(L1+L2) 

1020  L4=Ii2-L1 

1030  L5=0. 280675 1*L4 

1 040  T3=0 , 5*T,4*  ( PNC  ( L3+L5 )+ FMC  ( L3-L5 ) ) 

1050  T4='^XP(-Hrj]^^T3) 

1060  DEF  PND(D)=EXP(-D*3/4.1) 

1070  IX)=H3-H4 
1080  L1=H4 
1000  L2=n4+L0 
1100  L3=0.5*(L1+L2) 

1110  L4=L2-T,1 

1120  L5=0. 2886751 *L4 

1130  T5=0 . 5*L4*  ( PMD(  L3+L5  )+FNI)(  L3-L5 ) ) 

1140  T6=EXP(-0.128*T5) 

1150  T[J,2]=T4*T6 

1160  GOTO  1250 

1170  L0=:H3 

1180  L1=H0 

1190  L2=L0 

1200  L3=0.5*(L1+E2) 

1210  L4=L2-L1 

1220  L5=0. 2886751*^4 

1 230  T7=0. 5*E4* ( PND( L3+L5 )+Pra ( L3-L5 ) ) 

1240  Trj,2]=E]CP(-H[j3*T7) 

1250  IF  0=1  THEN  1270 

1260  PRINT  '»  HAZE/POG 

1270  REM  CALOirLATE  TRANSMITTANCE  OWING  TO  ATTENUATION  BY  RAIN. 

1280  IF  P=1  THEN  1310 

1290  Trj,3]=1 

1300  GOTO  1330 

1310  IF  V0>20  TITEN  1290 

1320  T[J,3]=EXr(-H3*Rrj] ) 

1330  IF  0=1  THEN  1350' 

1340  PRINT  " RAIN 

1350  REM  CALCULATE  TRAl^SMITTANCE  OV/ING  TO  ATTENUATION  BY  SMOKE. 

1360  T[j,4l=0.02/(Trj,l]*T[j-,2]*T[j,3]) 

1370  IF  t[j,4]  <=  1 THEN  1390 

1380  T[J,4]=1 

1390  IP  D=1  THEN  1420 

1400  PRINT  " SMOKE 

1410  PRINT 

1420  REM  CALCULATE  SMOKE  CONCENTRATION, 

1430  FOR  K=1  TO  4 

1440  REAL  D[K,1],D[K,2],D[K,3] 

1450  NEXT  K 

1460  IP  T[J,4]#1  then  1510 

1470  FOR  1=1  TO  4 

1480  C[J,I]=0 

1490  NEXT  I 

1500  GOTO  1570 


= ";T[J,2] 


'•;T[J,3] 


";T[J,4] 
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V 


y'itri 


1510  T8=L0G(T[J,4]) 

1520  T9=T8*T8 
1550  FOR  K®1  TO  4 

1 540  C[ J,K] =D[K, 1 ]+D[K, 2]*T8+D[K, 3]*T9 

1550  NEXT  K 

1560  IP  D=1  THEN  1610 

1570  PRINT  »'  SMOKE  CONCENTRATION: 

1580  PRINT  " 

1590  PRINT  " 

1600  PRINT  '• 

1610  NEXT  J . 

1620  IF  T>=^  THEN  1650 
1630  PRINT 
1640  PRINT 

1650  DISP  "DONE  - LINK  2" 

1660  DATA  0.118,26.7 
1670  DATA  0.22,13 
1680  DATA  0.14,5.3 
1690  DATA  0.55,5.1 
1700  -DATA  0,5 

1710  DATA  " 0.55  MICROMETERS:" 

1720  DATA  0.0093,-0.3428,-0.0009 
1730  DATA  0.0119,-0.2747,-0.0013 
1740  DATA  0.0142,-0.111,0.00004 
1750  DATA  0.0055,-0.1541,-0.0004 
1760  DATA  " 1.06  MICROMETERS:" 

1770  DATA  0.0093,-0.3428,-0.0009 
1780  DATA  0.0119,-0.2747,-0.0013 
1790  DATA  0.0142,-0.111,0.00004 
1800  DATA  0.0055,-0.1541,-0.0004 
1810  DATA  " 2.30  MICRCXIETERS ; " 

1820  DATA  0.0093,-0.3428,-0.0009 
1830  DATA  0.0119,-0.2747,-0.0013 
1840  DATA  0.0142,-0.111,0.00004 
1850  DATA  0.0055,-0.1541,-0.0004 
1860  DATA  " 3.80  MICROIffiTERS: " 

1870  DATA  0.0003,-0.3428,-0.0009 
1880  DATA  0.0119,-0.2747,-0,0013 
1890  DATA  0.0142,-0.111,0.00004 
1900  DATA  0.0055,-0.1541,-0.0004 
1910  DATA  " 10.6  MICRCMETERS : " 

1920  DATA  0.0093,-0.3428,-0,0009 
1930  DATA  0.0119,-0.2747,-0.0013 
1940  DATA  0.0142,-0.111,0.00004 
1950  DATA  0.0055,-0.1541,-0.0004 
I960  END 


FOG  OIL  - GM/SQ  M = 
HC  - GM/SQ  M = 
FS  - GM/SQ  M = 
V/P  - GM/SQ  M = 


";C 

";C 

";C 

";C 


J,1. 

J.2. 

J.5. 

[J,4] 
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COM  D,C0,C1,V0,P,T0,T1,D0,S0,Y,P0,R0,T[5,4],C[5,4],X0,Q[5*4] 

DIFFUSION  AND  SMOKE  SOURCE  STRENGTH  CALCULATIONS. 


10 

20  REM  m^IK:  ATMOSPHERIC 
30  Dm  S[6,4],A[6] 


40  DIM  A^[21 
50  FIXED  2 
60  IF  D=1  THEN  90 

70  PRINT  " ATMOSPHERIC  DIFFUSION  AND  SMOKE  SOURCE  STRENGTH  CALCULATIONS:" 
80  PRINT 

90  DISP  "TOTAL  DISTANCE  TO  BE  SMOKED  - M"; 

100  INPUT  XO 

110  DISP  "RELEASE  HEIGHT  OP  SMOKE  SOURCE"; 

120  INPUT  HO 

130  DISP  "MEAN  HEIGHT  OP  TARGET  - METERS"! 

140  riPUT  ZO 

150  DISP  "DIRECTION  OF  LINE  OP  SIGHT-DEG"; 

160  INPUT  AO 
170  IF  D=1  THEN  220 

180  PRINT  " TOTAL  DISTANCE  TO  BE  SMOKED  - METERS  = ";X0 

190  PRINT  " REIiEASE  HEIGHT  OF  SMOKE  SOURCE  ('\GL)  - METERS  = ";H0 

200  PRINT  " MEi\N  HIMGHT  OP  TARGET  - lUDTERS  = ";Z0 

210  PRINT  " DIRECTION  OP  LINE  OP  SIGHT  TO  TARGET  - DEG  = ";A0 

220  REM  DIPPUSION  CALCUT.ATIONS  POR  CONTII'IUOUS  SOURCE. 

230  POR  1=1  TO  6 
240  READ  A[I] 

250  r^XT  I 
260  POR  1=1  TO  6 
270  FOR  J=1  TO  4 
280  READ  STljJl 


290  NEXT  J 
300  JIEXT  I 

310  A1=-1.24+1 .19*LGT(Y) 

320  Z=10'‘A1 
330  D1=L0G(Z} 

340  B2=LOG(z)^2 
350  B3=LOG(z)'‘3 
360  B4=LOG(z)“4 
370  B5=L0G(Z)'‘5 

380  B6=0. 444685869+0. 294049265*B1-0. 23721 3914* B2 
390  B7=0. 1 55549504*B3-0. 03201 5723*B4+2. 1 51 68E-03*B5 
400  D1=B6+B7 
410  D1=EXP(D1  ) 

420  B6»-1 . 298283909-1 .006186784*B1+1 . 485094886*B2 

430  B7=-0.774136725*B3+0.156559355*B4-0.010823351*B5 

440  D2=B6+B7 

450  D2=EXP(D2)-0.225 

460  IF  Z>9. 999999  THEN  500 

470  B6=  5 . 77267E-04+2 . 31 943E-05*B1 +3 . 71 041 E-05*B2 
480  B7=-8. 40602E-06*B3+1 . 3421 E-07*B4+2. 551 31 E-08*B5 
490  GOTO  560 
500  IP  Z>40  THEN  540 


60 


51 0 B6=-1 1 . 56134901 +2 . 1 4824281 4*B1 -0 . 1 5621 0817*B2 
520  B7=7.03582E-03*B3-1.47353E-04*B4+1.18256E-06*B5 
530  (WTO  560 

540  B6=1 108. 366588-103. 5495836*B1+2.424499256*B2 
550  B7=-0. 01  -i584773*B3+4 . 3451 7S-05*B4-4. 695 56E-0S*B5 
560  I)3^B6+B7 

570  B6=0. 500775609+1 . C9261 4788*B1-1 . 573065836*B2 

530  B7=0 . 724276  579*B3-0 . 1 40820904*B4+9 .61621 E-03*B5 

590  D4=B6+B7 

60C  D4=EXP(D4)-1.2 

610  IE  Z>10  THEN  640 

620  B1  =L0(J(I)1  *X0*D2*1  / ( 1 +]:)3*XO‘‘D4  ) ) 

630  GOTO  650 

640  B1  =L0G(D1  *X0'‘r)2*(  1 +1  / (D3*X0"D4 ) ) ) 

650  B2=S[P0,1]*X0'‘S[P0,2V(1+S[P0,3J*X0*B[P0,4]  ) 

660  32=31 *B2 
670  S1=.4[PO]*XOW.9 
680  IP  3=1  THEN  740 
690  PRINT 

700  PRINT  " CONTINUOUS  SOURCE;" 

710  PRINT 

720  PRINT  " SIGMA  Y - METERS  = ";S1 

730  PRINT  " SIGM4  Z - METERS  = ";S2 

740  REM  SMOKE  SOURCE  STRENGTH  CALCULATIONS  FOR  CONTINUOUS  SOURCE. 
750  ,\2=ABS(AO-30)*(PI/180) 

760  IE  SOW  TH?:N  780 

770  S0=1 

780  S3=S0*0.515 

790  ^W=S2*S3+3QR(PI  )/SOR(  2 )*EXP(-0. 5*(  ( ZO-HO  )/S2  )'‘2 ) 

800  R2=S0R ( 3 . 3 1 24/ ( 6 . 76*S IN( A2 )*SIN( A2  UO . 49*C0S ( A2 )*COS ( A2 ) ) ) 
810  Y2=0.0337+(0.0369*R0)-(7E-04*R0*R0)+(6.1 1E-06*R0*R0*R0) 

820  Y3=1 . 3775+ ( 0 . 09868*R0 )- C 1 . 8E-03*RO*R0 )+ ( 1 . 56E-05*R0*R0*R0 ) 

830  Q0=r'0*R2 

840  Q1=00/Y2 

850  Q2=Q0/Y3 

860  FOR  1=1  TO  5 

870  ori,ll=C[l,l]*QO 

880  Q[l,2‘ =C ■l,2]*Q1 

890  Q[i,3]=C[I, 31*02 

900  BFJJ)  AS 

910  IF  3=1  THEN  900 

920  PRINT 

930  PRINT  AS 

940  PRINT 

950  PRINT  " SOURCE  STRENGTH:  FOG  OIL  - GMS/SEC  = ";Qri,1 

960  PRINT  " HC  - GMS/SEC  = ";Q  1,2' 

970  PRINT  " PS  - GMS/SEC  = ";Q[I,3 

980  NEXT  I 

990  REM  DIFFUSION  CALCULATIONS  FOR  QUASI- INSTANTANEOUS  SOURCE. 
1000  RESTORE  1350 
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1010 

1020 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1100 

1110 

1120 

1130 

1140 

1150 

1160 

1170 

1180 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 


31=2/3*31 
32=2/3*32 
IP  D=1  THEN  1090 
PRINT 

PRINT  " OUA3I-IN3TANTANE0U3  SOURCE:" 

PRINT 

PRINT  " 3IGMA  Y - METERS  = ";S1 

PRINT  " SIffl'IA  Z - METERS  = ";S2 

REM  SMOKE  SOURCE  STRENGTH  CALCULATIONS  FOR  QUASI- INSTANT AlffiOUS  SOURCE. 
00=( SI  *S2*PI )/EXP( - ( ZO*ZO-HO*HO )/ ( 2*S2*S2 ) ) 

Y4=3. 2469+ (0.0774*R0 )- ( 1 . 6E-03*R0*R0 )+ ( 1 . 73E-05*R0*R0*R0 ) 

00=(Qo/Y4)*10 

FOR  1=1  TO  5 

Q[l,4]=C[l,4l*Q0 

READ  AS 

IF  D=1  THEN  1210 
PRINT 
PRINT  AS 
PRIITT 

PRINT  " SOURCE  STRENGTH:  \fP  - GI^S/SEC  = ";Q[I,4] 

NEXT  I 


IP  D=1  THEN  1250 

PRniT 

PRHIT 

DISP  "DONE  - LINK  3" 

REM  DATA  USED  TO  CALCULATE  SIGI-IA  Y FOR  CONTINUOUS  SOURCE. 
DATA  0.4,0.32,0.22,0.144,0.102,0.076 

REM  DATA  USED  TO  CALCULATE  SI®IA  Z FOR  CONTINUOUS  SOURCE. 

DATA  0.112,1.06,5.38E-04,0.815 

DATA  0.13,0.95,6,52E-04,0.75 

DATA  O. 112, 0.92, 9. 05E-04, 0.718 

DATA  0.098,0.889,1.3535-03,0.688 

DATA  0.0609,0.895,1. 96E-03 ,0.684 

DATA  0. 0638,0. 783, 1.36E-03, 0.672 


DATA 

DATA 

DATA 

DATA 

DATA 

END 


0.55  MICRCMETSRS: 

1.06  MICROMETERS; 
2.30  MICROI-IETERS; 
3.80  MICRCIAETERS; 

10.6  MICROIffiTERS; 


10  COM  D,CO,C1,VO,P,TO,T1,DO,SO,Y,PO,RO,T[5,4],C[5,4],XO,Qr5,4] 

20  REM  ir./IK:  MUNITION  EXPENDITURES. 

30  DIM  Gr2,2],H[2,2],W[2,2],A[2,2],B[2,2],D[2,6],E[2,6],P[2,2],R[2,2] 
40  DIM  a:!5[21J 
50  EIXED  2 
70  PRINT  " 

80  PRIOT 


MUNITION  EXPENDITURES:" 


290 

300 

310 

320 

330 


480 

490 

500 

510 


100 

READ 

A 

1,1 

,A 

1,2 

,A 

2,1 

,A 

2,2. 

110 

READ 

B 

1,1 

,B 

1,2 

,B 

2,1 

,B 

.2,2. 

120 

READ 

D 

1,1 

,D 

1,2 

,D 

.1,3 

1,4. 

1,5 

,D 

1,6. 

130 

READ 

D 

2,1 

,I> 

2,2 

,0 

2,5 

,D 

.2,4] 

,D 

.2,6 

140 

READ 

E 

1,1 

,E 

.1,2 

.1,3 

,E 

,E 

.1,5 

1,61 

150 

READ 

jL$ 

[2,1] 

2,2] 

,15 

[2,3 

,E 

2,4] 

,35 

[2,5l 

,33 

2,6] 

160  DISP  "TME  SMOKE  RE(1UIRED  - MINUTES"; 
170  INPUT  T2 
180  FOR  K=1  TO  5 
190  IF  K=1  THEN  210 
200  IP  D=1  THEN  830 
210  IF  T[K,4l=1  THTN  760 
220 


230 

Ct 

1,1 

3^,2 

1,1 

240 

Ct 

11,2. 

=0 

.3-.  4 

/H 

1,2. 

250 

a 

.2,1. 

K,2 

2,1. 

260 

Q 

[2,2] 

=Q 

[3^4] 

/H 

2,2] 

270  FOR  1=1  TO  2 
280  FOR  J=1  TO  2 
GO=TKT(G[T, J] ) 
G1=0[T,JUG0 
IP  01=0 


330 


THEN 

or  I, J] =00+1 
IWXT  J 
340  NEXT  I 

350  REM  CALCULATE  TOTAL  TIME  FOR  REPLENISH-IENT. 

360  FOR  1=1  TO  2 
370  FOR  J=1  TO  2 
380  V/ri,j]=T2+A[l,j]-B[l,J] 

NEXT  J 
NEXT  I 

REM  CALCULATE  RATE  OF  FIPJ:. 

420  FOR  J=1  TO  2 
430  F[l,Jl=V.'ri,jl*D[j,IX)l+1 
440  Fr2,Jj=\7[2, J]*E[J,P0]+1 
450  I'TEXT  J 

460  IISI4  CALCULATE  TOTAL  NUMBER  OF  ROUNDS  REQUIRED, 
470  FOR  1=1  TO  2 
FOR  J=1  TO  2 


390 

400 

410 


R[l,J]=Gri,Jl 

R1=INT(Rri,Jj 

R2=R[I,J]-R1 


*P[I,J] 
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i 


520  I?  R2=0  THEN  540 
530  R[I,J]=H1+1 
540  ICXT  J 
550  NEXT  I 
560  RE.M)  A:; 

570  IP  n=1  THJCN  6-y) 

580  PRINT  A?.; 

590  PRINT 

600  PRINT  " VOTiUI'IE  OP  PIRJC  - HC  SKOKES GREEN” 

610  PRINT  ” SCREEN  LENGTJI  - I-!ETERS  = ”;X0 

620  PRINT  ” SCREEN  Dmi.’TION  - XrRJTES  = ";T2 

630  PRINT 

640  PRINT  " GUNS  ROUNDS/lIIN  TOTAL  ROUNDS” 

650  PRINT  ” 105  ";G[l,l];''  ”;D[l,POl;”  ";H[1,i1 

660  PRINT  " 153  ";G[2,l];”  ”;E[i,POJ;”  ";R[2,i] 

670  PRINT 

630  PRINT  " VOLIRffi  OF  PIRl^  - VI?  SMOKESCREEN” 

690  PRINT  " SCREEN  LENGTH  - !'!ET>DR3  = ”;X0 

700  PRINT  ” SCREEN  DURATION  - MIIUTES  = ”;T2 

710  PRINT 

720  PRINT  " GUNS  ROUNDS/MIN  TOTAL  ROUNDS” 

730  PRINT  ” 105  ";G[l,2l;”  ”;Dr2,P0];” 

740  PRINT  " 155  ";Gr2,2];”  ";Er2,P0l;"  ";R[2,2] 

750  GOTO  010 

760  READ  AO 

770  IP  D=1  THi^N  800 

780  PRINT  AO 

790  PRINT 

800  PRINT  ” SMOKE  NOT  REQUIRED  DUE  TO  ATMOSPHERIC  OOTIDITIONS.  ” 

310  PRINT 
820  NJJXT  K 
830  PRINT 
840  PRINT 
350  DT3P  "DONE” 

860  REM  UNIT(PER  GUN)  SOURCE  STRENGTHS. 

870  DATA  18.9,1737.3,46.0,7076.2 
880  REM  TIME  FOR  BUILDUP. 

890  DATA  1,0. 5, 1,0. 5 
900  REM  AVERAGE  BURN  TIME. 

910  DATA  3,0.0167,4,0.0167 

920  REM  PuATE  OP  FIRE  VS.  STABILITY  CATEGORY, 

930  DATA  6, 4, 3, 2, 1,1 

940  DATA  0,0, 6, 4, 1.5,1 

950  DATA  3,2,1.5,1,0.5,0.333 

960  DATA  0,0,3,2,0.5,0.333 

970  DATA  " 0.55  MICROIffiTERS : " 

980  DATA  ” 1.06  MICROMETERS:" 

990  DATA  ” 2.30  MICROMETERS:” 

ioOO  DATA  ” 3o80  MICRC^IETERS:  ” 

1010  DATA  ” 10.6  MICROffiTERS:” 

1020  END 


HPL 

GLOSSARY  OF  MNEMONICS 


A 

B 

C 

D 

E 

F 

G 

H 

P 

Q 

R 

T 

Y 

Z 

A(7,9) 

B(5,4) 

C(5.4) 


D(5) 

E(5) 

F(5) 

G(5) 

H(3,4,5) 


1(5,4) 

J(6) 

K(6,4) 

M(4) 
P(2,2) 
0(2,2) 
R(2,2) 
S(2,2) 
T(2,2) 
U(2,2) 
V 2,6) 

WU,6) 

X(2,2) 

A$(3) 

B$(6) 

C$(80) 

D$(3) 

E$ 


Ceiling  - hundreds  of  feet 
Cloud  cover  - per  cent 
Visibility  - miles 
Temperature  - degrees  F 
Dew  Point  - degrees  F 
Wind  dirction  - tens  of  degrees 
Wind  Speed  - knots 
Atmospheric  stability  category 
Mixing  depth  height  - meters 
Relative  humidity  - percent 
Total  distance  to  be  smoked  - meters 
Time  smoke  required  - minutes 
Average  roughness  element  - centimeters 
Roughness  length  - centimeters 

Table  of  stability  categories  (depending  upon  solar 
altitude  and  wind  speed) 

Table  of  transmittances  owing  to  water  vapor,  haze/fog. 

rain  and  smoke  for  0.55,  1.06,  2.3,  3.8,  and  10.6  micrometers 

Table  of  smoke  concentration  values  for  fog  oil,  HC,  FS  and 

WP  for  0.55,  1.06,  2.3,  3.8,  and  10.6  micrometers 

Error  function  absorption  coefficients 

Scale  height  for  Mie  scattering 

Haze  and  fog  attenuation  coefficients 

Rain  attenuation  coefficients 

Table  of  coefficients  used  to  calculate  smoke  concentrations 
using  the  calculated  transmittance  values  for  0.55,  1.06,  2.3, 
3.8,  and  10.6  micrometers 

Smoke  source  strength  values  for  fog  oil,  HC,  FS,  and  WP 

for  0.55,  1.06,  2.3,  3.8,  and  10.6  micrometers 

Coefficients  to  compute  ay  - continuous  source 

Coefficients  of  the  roughness  correction  factor  used  in 

calculating  az  for  the  various  roughness  lengths 

Yield  factors  for  fog  oil,  HC,  FS,  and  WP 

Total  number  of  rounds  required  to  maintain  smoke  screen 

Number  of  guns 

Number  of  rounds  per  gun 

Unit  (per  gun)  source  strength 

Smoke  build-up  time 

Munition  average  burn  time 

Rate  of  fire  vs  stability  category  for  105  Howitzer 

Rate  of  fire  vs  stability  category  for  155  Howitzer 

Total  time  for  munition  replenishment 

Met  observation  station  identifier 

Stability  category  indicator 

Wavelength  indicator 

Precipitation  indicator 

Demo  indicator 
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1:  dim  A [7 , 9 ] ,fl [5,4 1 ,C(5, 4) ,D [5 ] , E (5 ) , F [5  I ,G(51 , H( 3,4 , 5] , I [ 5, 4) 
2:  dim  J16) ,K16,41  ,M[41 ,P12,21 ,Q12,21  ,R[2,21 
3:  dim  S [2, 21 ,T [ 2, 21 , U[ 2, 2] , V[2,6 J , W(2,6 ] , X [2, 2] 

4:  dim  A$ [3 ] ,B$  [6]  ,C$ [ 35) ,D$ [31 , E$ (3 J 
5;  asgn  "KDATAl",! 

6:  files  KOATAl 

7:  sread  i, A [*  1 , D [* 1 , E ( *1  , H[ *1  , J [ *1 , K [*  1 
8;  sread  1 ,3 [* 1 , T [* 1 ,U [ * 1 , V [*1  , W [* 1 
9;  sread  1,8$, C$ 

10;  fmt  2/;wrt  701 

11;  fmt  "KWIK  SMOKE  PROGRAM" ;wrt  701 
12:  fmt  /;wrt  701 

13:  ent  "IS  THIS  A DEMO?  YES  OR  i40",E$ 

14:  ent  "MET  SITE  ID",A$ 

15;  ent  "LATITUDE  OF  MET  SITE  - DEG",r2 
16:  ent  "LONGITUDE  OF  MET  SITE  - DE3",r3 
17;  ent  "ALTITUDE  OF  MET  S ITE-KILOMETERS" , r 4 
18;  ent  "JULIAN  DATE  OF  MET  OBSERVATION" , r 5 
19:  ent  "ZULU  TIME  OF  MET  OBSERVATION-HR" , r 6 


20: 

fmt 

If 

MET  SITE:";wrt 

701 

21; 

fmt 

n 

";wrt  701 

22: 

fmt 

n 

ID 

= ",2x,c3;wrt  701, A$ 

23: 

fmt 

It 

LATITUDE 

- DEG 

= ",f6.2;wrt  701, r2 

24: 

fmt 

It 

LONGITUDE 

- DEG 

= ",f6.2;wrt  701, r3 

25; 

fmt 

II 

ALTITUDE 

- KM 

= " ,2x,f4.2;wrt  701, r4 

26: 

fmt 

II 

";wrt  701 

27; 

fmt 

II 

JULIAN  DATE 

- DAY 

= ",f3.0;wrt  701, r5 

28: 

fmt 

11 

ZULU  TIME 

- HOUR 

= ",lx,f2.0;wrt  701, r6 

29: 

fmt 

/; 

wrt  701 

30: 

ent 

"CEILING  - HUNDREDS 

OF  FEET 

" ,A 

31: 

A*100* 

.3043-A 

32:  ent  "CLOUD  COVER  - PERCENT" ,B 
33:  ent  "VISIBILITY  - MILES" ,C 
34:  C*1.61*C 

35;  ent  "PRECIPITATION  - YES  OR  NO",D$ 

36;  ent  "TEMPERATURE  - DEG  F",D 
37:  5/9*(D-32)-D 
38:  ent  "DEN  POIiNT  - DEG  F",E 
39;  5/9* (E-32)*E 

40:  ent  "WIND  DIRECTION  - TENS  OF  DEGS",F 
41:  F*10*F 

42:  ent  "WIND  SPEED  - KNOTS", 3 

43;  ent  "AVE  ROUGHNESS  ELEMENT  - CM",Y 


44: 

fmt 

II 

METEOROLOGICAL  INPUTS: 

";wrt  701 

45: 

fmt 

»l 

";wrt  701 

46: 

fmt 

If 

CEILING 

- METERS 

S 

",f8.2;wrt  701, A 

47: 

fmt 

II 

CLOUD  COVER 

- PERCENT 

a 

",f8.2»wrt  701, B 

43: 

fmt 

It 

VISIBILITY 

- KILOMETERS 

a 

",f8.2;wrt  701, C 

49: 

fmt 

M 

PRECIPITATION 

« 

*,4x,c3;wrt  701, D$ 

50; 

fmt 

It 

TEMPERATURE 

- DEG  C 

» 

",f3.2;wrt  701, D 
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51: 
52: 
53: 
54: 
55: 
56: 
57: 
53: 
59: 
60: 
51: 
62: 
63: 
64 : 
55: 


E.nt  " DEWPOIHT  - DEG  C 

fmt  " direction'!  - DEG 

fmt  " vifli^D  SPEED  - KNOTS 

fmt  " AVE  ROUGHNESS  ELEMENT  - CM 

fmt  /;wrt  701 
if  B#100;gto  "KIOOO" 
if  A>2133.6042 ;gto  "KIOOO" 

0*r0 

0*rl 

gto  "K1400" 

"KIOOO": 

fmt  "CALCULATE  ANGULAR  FRACTION  OF  A YEAR  FOR 
( r5-l)*360/365.242*r9 

fmt  "CALCULATE  SOLAR  DECLINAriON  ANGLE" 
279.9348+r9-rll 


» ",f8.2;wrt  701, E 
» ■,f8.2;wrt  701, F 
» ",f3.2;wrt  701, G 
- ",f8.2;wrt  701,y 


A GIVEN  JULIAN  DATE" 


66:  rll+1.914327*sin(r9)-.079525*co3(r9)*rll 

67:  rll  + .019938*sin(2*r  9)-.00162*cos(2*r9)Tll 

68:  23.4433-*rl2 

69:  sin { rl2) *sin ( r 11  )*r 13 

70:  a3n(rl3)*rl3 

71:  fmt  "CALCULATE  TIME  OF  MERIDIAN  PASSAGE  - TRUE  SOLAR  NOO^" 

72:  12  + . 12  357*sin(:x9)-.004289*cos(c9)  ♦rl4 

73:  rl4  + . 153809*3 in{2*c 9) +.060783*c»s  (2*c 9) ♦clM 

74:  fmt  "CALCULATE  SOLAR  HOUR  ANGLE" 

75:  15* ( r 6-rl4) -r 3*rl5 

76:  fBt  "CALCULATE  SOLAR  ALTITUDE" 

77:  sin(r2)*sin(rl3)+cos(r2) *co3(rl3) * cos ( r 15 ) ♦r 16 
78:  a3n(tl6)*rl6 

79:  fmt  "CALCUl,ATE  TIME  OF  SUNRISE  AND  SUNSET" 

80:  -1.76459*r4'.40795*rl7 

81:  (sin(rl7)-sin(r2)*sin(ri3)) / (cos (r  2) * cos( rl3) ) *r 18 

82:  acs(rl8)*rl8 

83:  rl8* (24/360)*rl8 

84:  r3/15ftl4-rl3*tl9 

85:  r3/15frl4  + rl8T20 

86:  if  r 20>24 ;r20-24*r20 

87:  fmt  "CALCULATE  INSOLATION  CLASS  NU4BER" 

38:  0*cl 

89:  if  rl6>60;4*rl;gto  "KllOO" 

90:  if  tl6>35;3Tl;qto  "KllOO" 

91:  if  t 16>15 ;2*rl ;gto  "KllOO" 

92:  if  rl6<=*0;qto  "K1300" 

9 3:  l*tl 
94.  "KllOO": 

95:  fmt  "CALCULATE  NET  RADIATION  INDEX  FOR  DAYTIME." 

96:  0t2 

97:  if  3<50;r l*r2;gto  •K1200" 

98:  if  A<2i33.6042;rl-2*r2;gto  "K1200" 

99:  if  A<4376.8096;rl-2*r3;ito  "K1200" 

100:  if  3*100;rl-l*r 2 
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101;  "K1200": 

102;  if  t2=0;rl>r2 
103;  if  r2<l;lT2 
104;  r2*r0;gto  "K1400" 

105;  "K1300"; 

108;  fmt  "CALCULATE  MET  RADIATION  INDEX  FOR  MIGHTTIME." 
107;  if  3<40;-2*r0;gtO  "K1400" 

108;  -l»r3 
109;  "K1400"; 


110; 

fmt  "CALCULATE  PA3QUILL  STABILITY  CATEGORIi." 

111; 

O- 

c4  ;0*r  5 

112; 

if 

r0=4  ;i-*'r4 

113; 

if 

r 0 *3 ; 2*r4 

114; 

if 

r0=2;3-r4 

115; 

if 

r 0=1; 4>r4 

116; 

if 

r0=0;5*r4 

117; 

if 

r 0 =-l ; 6*r 4 

118: 

if 

r0=-2;7*r4 

119: 

if 

G<2;1t5 

120: 

if 

G<4  ;2t5 

121: 

if 

G<4;2>r5;gto  " 

K1500'’ 

122: 

if 

G<6;3*r5;gto  " 

K1500" 

123: 

if 

G<7;4*r5;gto  " 

K1500" 

124: 

if 

G<8;5*r5;gto  " 

K1500" 

125; 

if 

G<10 ;6*r5 ;gto 

"K1500" 

126: 

if 

G<11; /♦r5 ;gto 

"K1500" 

127: 

if 

G<12 ;3*r5 ;gto 

"K1500" 

128: 

9* 

r5 

129;  "K1500"; 

130;  Alr4,r5]-H 

131;  fmt  "CALCULATE  MIXING  DEPTH  HEIGHT." 

132;  (6-U) "121* (D-E)/6+H*.087* (G+.5)/(12*8.237e-5*5.809)*P 
133;  fmt  "CALCULATE  RELATIVE  HUMIDITY." 

134;  if  D>0;gto  "XloOO" 

135;  9.5T0;265.5*rl 
136:  gto  "X1700" 

137;  "K1600": 

138;  7.5t0;237.3t1 
139:  "K1700"; 

140;  if  E>0;gto  "K1800" 

141:  9.5*r2;265.5*r3 
142;  gto  "K1900" 

143:  "K1800"; 

144:  7.5*t2;237.3*r3 
145:  "K1900": 

146;  6.11*10“(r0*D/(rl+D))T4 
147;  6.11*10'“(r2*E/(r3+E)  )t5 
148;  r5/r4*100*Q 

149;  fmt  " METEOROLOGICAL  CALCULATIONS ;"; wr t 701 

150;  fmt  " ";wrt  701 
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151;  fnit  ’’  PASQUILL  3rA3ILIT5r  CATEGORY  » ",5x,cl 

152:  wrt  701,B$[tI,Hl 

153;  fmt  “ RSLATIv/E  dUMIDIfY  = ",f6.2 

154;  wrt  701,0 

155:  if  £$  (1,3]  = "YES"  ;gto  "jC2000" 

156:  fmt  /;wrt  701 

157:  fmt  " AiMOSPHERIC  OPTICS  A«D  S.<OKE  CORCErJTRATION  CALCULATIORS : " 
153:  wrt  701 
159;  "K2000": 

160;  fmt  •'  ";wrt  701 
161:  ln(C)>r3 
162:  r0*r0*rl 
163:  rl*r0*r2 

164:  1.5551-.9311*r0-.0197*rl+.0041*r2*F(l] 

165;  axD(F(ll)*Fll] 

166:  1.5551-.-9311*rO-.0197*rl  + .004  1*r2*F[21 
167:  exo{P(2] )^F(21 

168:  1. 44 91-1. 0044 *r0-. 012* rl+. 0032* r2*F[ 31 
169;  exo(F[31 )>f (31 

170:  1.2i94-L.’04  36*rO  + .0099*rl-.0016*r2*F[41 
171:  exp ( P [4 1 ) *F(4 1 

172;  1.5176-1.714  7*  rO 0001*  r 1+.  0428*r  2*F( 5 ) 

173:  exp(F(51)*P[51 

174;  1. 33  06-.3825*  1 0- . 07  53* rl  + . 012 9*r2*a ( 1 ) 

175;  exD(3(l])*G(l] 

176:  1.4  00d-.9  36  5*r0-.  J14*rl4-2.  3e-3*r2*G(21 
177;  exp(0(2] )*G(2] 

173:  1.5497-.3696*r0-.1084*rl+.0231*r2*Gl31 
179;  exp(3[31 )*G(3] 

130;  1.5556-.9013*r0-.0773*rlf.0173*r2*3l41 
181;  exp{G(4] )*a(4J 

182:  1.5928-.9396*r0-.0627*rl+.0168*r2*G(5) 

183;  exp(G(5| )+G(5) 

134;  ent  "SLANT  RANGE  TO  TARGET  - METERS" ,r20 
135;  ent  “ANGLE  OF  SIGHT  TO  TARGET  - DEG",r6 
186:  if  E$[l,3i="YE3";gto  "K2100" 

187:  fmt  " SLANT  RANGE  TO  TARGET  - METERS  » ",f8.2 

138;  wrt  701, t20 

189:  fmt  " ANGLE  OF  SIGHT  TO  TARGET  - DEG  » ",3x,f5.2 

190;  wrt  701, r6 
191;  "K2100": 

192:  r20/1000*r20 
193:  if  r6<0;-r6*r6 
194;  sin(r6)*r6 
195;  0*r8 

196;  if  r6=*0;gto  "K2200'’ 

197;  l/r6*r8 
193:  "K2200"; 

159;  fmt  "CALCULATE  PRECI PITAiLE  WATER." 

200;  .4477+.0323*E+1.2e-3*E*2+1.84e-5*E*3*rll 
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201:  if  E$[l,3]  = “YES";gto  •j<2300'‘ 

202:  fmt  " ";wrt  701 

203:  fmt  *•  PSRCIPITABLE  WATER  - CM/KM  * ■,4x,f4.2 

204:  wrt  701, cll 
205;  "K2300"; 

206;  fmt  "CALCULATE  A40UNT  OF  WATER  VAPOR  IN  PATri." 

207;  r20*r0;0*rl;t0*r2;.5*(rl+r2)^r3 
208:  r2-rlT4  ;.2886751*r4*r5 
209:  .5*r4*  ( 'PNA  ' (r  3+r5)  +'FNA'  {r3-r5)  )t9 
210;  rll*r9*rl0 

211:  if  E$ (1,3 ]="YES" ;gto  "K2400" 

212;  fmt  “ AMOUNT  OF  WATER  VAPOR  IN  PATH  - CM  » ”,4x,f4.2 

213:  wrt  701, rlO 
214:  "K2400"; 

215;  fmt  "CALCULATE  TRANSMITTANCE  FOR  0.55,1.06,2.3,3.8,10.6  MICROMETERS." 

216;  1*M;17*N 

217;  for  1*1  to  5 

21d:  if  E$ (1 ,3 ]*" ^ES" ;gto  "K2500" 

219:  fmt  " ";wrt  701 

220;  fmt  4x,cl7;wrt  701,C$[M,N] 

221:  "R2500"; 

222;  N+17*N;M+17-M 

223:  fmt  "CALCULATE  TRANSMITTANCE  OWING  TO  ABSORPTION  BY  WATER  VAPOR." 

224;  if  I*5;exp(-.0681*rll)^B[I,ll;gto  "K2600" 

225;  D(Il*/(rl0*ff)/2-t0;0-rl;r0-r2 
226:  .5"(rl+r2)*r3 
227:  r2-rl-r4 
228:  .2836751"r4*r5 

229:  . 5*r4* ('FNB'(r3  + r5) +'FN3  ' (r3-r5) )*rl2 
230:  2/|jT*rl2*a(I,  II 
231:  1-B[I,11»B(I,1) 

232;  "K2600": 

233;  if  E$[l,31*"yES" ;gto  "K2700" 

234:  fmt  " ";wrt  701 

235;  fmt  " TRANSMITTANCE  OWING  TO  ATTENUATION  BT:  WATER  VAPOR  - ",f5.2 

236:  wrt  701,B[I,1] 

237:  "K2700": 

2 33:  fmt  "CALCULATE  TRAi4SMITTANCE  OWING  TO  ATTENUATION  BY  HAZE  AND  FOG." 

239:  if  D$(L,3)*"YES";1*B[I,21 ;gto  "K2900" 

240:  if  C>»E(I)jgto  "K2800" 

241;  r 8^r0 ; 0*r 1 ; r0*r 2 
242;  .5*(rl+r2)*r3 
243:  r2-rl*r4 
244:  .2886751*r4*r5 

245:  .5*r4*('FNC'(r3+r5)+'FNC' (r3-r5) )*rl3 

246:  exp (-F [ I] *rl3) ♦r 14 

247;  r20-r8*r0 ;r3*rl ;t8+r0»r 2 

248;  .5* (rl+r2) ♦r3 

249:  r2-tl-t4 

250:  .2886751*r4*r5 
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251;  .5*r4*  {'Ft<D' (t3+t5) +'FNO'(r3-r5)  )*rl5 

252;  exp(-.128*rl5)Tl6 

253;  rl4*rl5*B{  1,2]  ;gto  "iC2900" 

254;  “K2300"; 

255;  c20*r0;0Tl;r0*r2;  .5*  (rl+c2)*r3 
256;  r2-rl*r4; .2886751* r4*r5 
257:  .5*r4*('FL'JD'(r3fc5)+'FMD'(r3-r5)  )>rl7 
258:  exp(-FlIJ*rl7)*3[I,2] 

259;  "K290a"; 

260;  if  e$ll,3|  ='"YES";gto  "K3000* 

261;  f^t  " HAZE/FOa 

262:  wet  701,3(1,21 
263;  "KOOOO": 

264:  f(tit  "CALCULATE  TRANSMITTANCE  OWING  TO  ATTENUATION  BY  RAIN." 
255:  if  D? (1,21  ="NO";l*B(I, 3J ;gto  "K3100" 

256;  if  C>20;l*a(I ,31 ;gto  "K3100" 

257;  exp(-r20*3(Il )*a(I,31 
268:  "i<3100": 

269;  if  E$ll,3|  ="YE3" ;gto  "K3200" 

27  0;  f:nt  " RAIN 

271;  wrt  701,3(1,31 
272:  "K3200“: 

27  3 : fmt  "CALCULATE  TRANSMITTANCE  OWING  TO  ATTENUATION  BY  SMOeCE." 
274:  .02/ (3 [I, 11 *3(1,21 *3 (1,31) *3 (I, 4] 

275:  if  3(1,41 >l;l*all,41 

276:  if  E$(l,31  ="YES"  ;gto  "K3300'* 

27  7;  fmt  " SMOKE 

278;  wet  701,a(I,4| 

279;  "K3300": 

230:  fmt  "CALCULATE  SMOKE  CONCENTRATION. " 

281:  if  8(1,4] #l;gto  "K3400" 

282:  foe  0*1  to  4;0*C(I,J1 
283;  next  J;gto  "K3500" 

284;  "K3400": 

285:  ln(B(I,41 ) ♦cia 
286;  cl8*el8*cl9 


287;  for  K=1  to  4 

288:  ri(l,K,I)  +H(2,K,  I]  * r 18+11  ( 3 , K, 1 1 * rl9*C(  I , K ] 
289:  next  K 
290:  "K3500"; 


291; 

if  E 

;$(i 

,31*’ 

•YES" ;gto  "K3600" 

292: 

fmt 

n *1 

jwr  t 

701 

293: 

fmt 

N 

SMOKE  CONCENTRATION: 

FOG  OIL 

- GM/SQ 

M 

294: 

wet 

701 

,C(I, 

,11 

295: 

fmt 

M 

HC 

- GM/SQ 

M 

296: 

wrt 

701 

,C(I, 

,21 

297; 

fmt 

II 

FS 

- GM/SQ 

M 

298; 

wrt 

701 

,C(I. 

,3] 

299: 

fmt 

M 

WP 

- GM/SQ 

M 

300; 

wrt 

701 

,C(I, 

,41 
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301:  "K3600": 

302:  next  I 

303:  iE  E$  (1 ,3  ] = " YES"  ;gto  "rC4000" 

304:  fmt  /;wtt  701 

305:  fmt  " ATMOSPHERIC  DIFFUSION  SOURCE  STRENGTH  CALCULATIONS wr t 701 
306:  fmt  " ";wrt  701 
307:  “K4000": 

308:  ent  "TOTAL  DISTANCE  TO  3E  SMOKED  - M",R 
309:  ent  "RELEASE  HEIGHT  OF  SMOKE  SOURCE", rO 
310:  ent  " 4EAN  HEIGHT  OF  TARGET  - METERS", rl 
311:  ent  "DIRECTION  OF  LINE  OF  SIGHT  - DEG",t2 
312:  if  E$ (1,3]="YES" ;gto  "K4050" 

313:  fmt  " TOTAL  DISTANCE  TO  BE  SMOKED  - METERS  = ",f8.2 

314:  wrt  701 , R 

315:  fmt  " RELEASE  HEIGHT  OF  SMOKE  SOURCE  (AGL)  - METERS  = ",2x,f6.2 

316:  wrt  701, rO 

317:  fmt  " MEAN  HEIGHT  OF  TARGET  - METERS  * ",2x,f6.2 

318:  wrt  701, rl 

319:  fmt  " DIRECTION  OF  LINS  OF  SIGHT  TO  TARGET  - DEG  =*  ",2x,f6.2 

320:  wrt  701, r2 
321:  "K4050": 

322:  fmt  "DIFFUSION  CALCULATIONS  FOR  CONTINUOUS  SOURCE." 

323:  -1.24+1. 19*log(Y)-r3 
324:  lO^rS-Z 

325:  ln(Z)-r8;ln(Z)"2*r9; ln(Z) "3*rl0 
326:  ln(Z)'4-rll;ln(Z)*5*rl2 

327:  .4 44635869+. 29404 926 5*r 3-. 2 37213914* r9*rl3 
328:  . 155349504*rl0-.032015;23*rll+2.15163e-3*rl2*rl4 
329:  rl3+rl4*r4 
330:  exp(r4)*r4 

331:  - 1.29828 3 909-1. 006 186784*r 8+1. 48 5094886*r9^r 13 

332:  -. 7 74 13672 5* rl0+. 1565593 55* rll-. 010 823351* rl2*r 14 

333:  rl3+rl4-*-r5 

334:  exp  (r5)-.225>r5 

335:  if  Z>9.999999;gto  "K4100" 

336:  5.77267e-4+2.31943e-5*r8+3.71041e-5*r9*rl3 
337:  -3.40602e-6*rl0+1.3421e-7*rll+2.5513ie-8*rl2*rl4 
338:  gto  "K4175" 

339:  "K4100": 

340:  if  Z>40;gto  "K4150" 

341:  -11. 56134901+2. 143242814*r8-.156210817*r9*rl3 
342:  7.03582e-3*tl0-1.47353e-4*rll+1.18256e-6*rl2-rl4 
343:  gto  "K4175" 

344:  "K4150": 

345:  1108.366588-103.5495836*r8+2.424499256*r9*rl3 
346;  -.014 584773*r 10+4. 3451 7e-5*r 11-4. 69 556e-8*rl2*r 14 
347;  "K4175"; 

343:  rl3+rl4-r6 

349;  . 500775609+1. 0926 14788* r8-l. 573065836* r9>r 13 
j 350:  .724276579*rl0-.140820904*rll+9.61621e^3*rl2*rl4 
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1- 


' , f3 .2; wrt 
' , £3 . 2 ; wet 


701, rll 
701, clO 


351:  rl3+rl4-r7 
352:  exp(r7)-1.2T7 
353:  if  Z>10;gto  "K4200" 

354:  ln(r4*R~r5*l/(l+r6*R"r7)  )t3 
355:  gto  "K4225" 

356:  "K4200": 

357:  ln{r4*R“r5*  (l+l/(  r6*R'‘r  7)  ) )*r8 
358:  ’•K4225": 

359:  K(8,l]*R"K(H,21/(l+KtH,31*R'‘K(H,41  )^r9 
360:  r8*r9*cl0 
361:  J[H)*R''.9*rll 
362:  if  E$  [1,3]  ="  YES"  ;gto  ■•K4250" 

363:  fmt  " ";wtt  701 

364:  fmt  " CONXINUOJS  SOURCE:"; wrt  701 
365:  fmt  " ";wrt  701 
366:  fmt  " SIG-IA  Y - METERS  = 

367:  fmt  " SIGMA  Z - METERS  = 

368:  "K4250": 

389:  fmt  "SMOKE  SOURCE  STRENjrU  CALCULATIONS  FOR  CONTINUOUS  SOURCE.' 

370:  abs ( r2-F) ♦rlZ 

371:  if  G=0;1-G 

372:  G*.515*rl3 

373:  1>M;17*N 

374:  cl0*rl3*/n//2*exp(-.5*  ( (rl-rO)  /tlO)  ''2)t14 

375:  /(3.3124/(6.76*sin(rl2) *sin (c 12) +. 49*cos( cl2) *cos ( cl2) ) )-rl5 

3 76:  .9337  + .0  369*Q-7e-4*Q*Q+6.'lle-6*Q*Q*Q*M[21 

377:  1. 3775  + . 09368*0-1. 8e-3*Q*Q+l. 56e-5*Q*Q*Q*M [3  I 

378:  rl4*rl5Tl4 

379:  rl4/M[2]Tl6 

330:  rl4/M[3j*rl7 

381:  for  K=1  to  5 

332:  C(K,l]*rl4*I [K,l] 

333:  C[K,21*rl6^I [K,21 
384:  C[K,31*rl7*I (K,31 
385:  if  E$[l,3]="YeS";gto  "K4300" 

336:  fmt  " ";wrt  701 

387:  fmt  4x,cl7;wrt  701,C$[M,N] 

388:  "K4300": 

389:  N + 17-N;:’1  + 17*M 

390:  if  E$ [1,3] ="YES" ;gto  "K4400" 

391:  fmt  " ";wrt  701 
392:  fmt  " SOURCE  STRENGTH: 

393:  wrt  701,I[K,1] 

394:  fmt  " HC 

395:  wrt  701, UK, 2] 

396:  fmt  " FS 

397:  wrt  701,I[K,3] 

398:  "K4400": 

399:  next  K 

400:  fmt  "DIFFUSION  CALCULATIONS  FOR  QUASI- INSTANTANEOUS  SOURCE.' 
*14590 


FOG  OIL  - GMS/SEC 

- GMS/SEC 

- GMS/SEC 


',£8.2 

',£3.2 

',£3.2 
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401:  2/3*rllTll 
402:  2/3*rl0-rl0 

403:  if  fc:$(l,3]="YES";gto  "K4500" 

404:  frat  " ";wrt  701 

405:  fmt  " QUASI- INSTANTANEOUS  SOURCE: ”;wrt  701 
406:  fmt  " ";wrt  701 

407:  fmt  " SIC^  Y - METERS  = ",f8.2;wrt  701, rll 

403:  fmt  " SIGMA  Z - METERS  = ",f8.2;wrt  701, rlO 

409:  ''K4500": 

410:  fmt  "SMOKE  SOURCE  STRENGTH  CALCULATIONS  FOR  QUASI-INSTANTANEOUS  SOURCE." 
411:  1*M;17*N 

412:  rll*rl0*Ti/exp(-  (r  1* r 1-r 0* rO ) / ( 2*cl0* rlO ) ) *rl4 
413:  3.2469+.0774*Q-1.6e-3*C*O+1.73e-5*Q*Q*Q-M[4] 

414:  cl4/M(41*10-rl4 
415:  for  J=1  to  5 
416:  C[J ,4]*rl4-I [ J,4] 

417:  if  E$ [1 ,3]="YES" ;gto  "K4600" 

418:  fmt  " ";wrt  701 

419:  fmt  4x,cl7;wrt  701,C$[M,N] 

420:  "K4600": 

421:  N+17+N ;M+17-M 

422:  if  E$[l,3]="YE3";gto  "K4700" 

423:  fmt  " ";wrt  701 

424:  fmt  " SOURCE  STRENGTH:  WP  - GMS/SEC  = ",f3.2 

425:  wet  701,I[J,41 
426:  "K4700": 

427:  next  J 
428:  fmt  /;wrt  701 

429:  fmt  " MUNITION  EXPENDITURES :" ;wr t 701 
430:  fmt  " ";wrt  701 

431:  ent  "TIME  SMOKE  REQUIRED  - MINUTES" ,T 

432:  1*I;17*J 

433:  for  K=1  to  5 

434:  if  K=l;gto  "K5000" 

435:  if  E$ll,31 ="YES" ;gto  "K5400" 

436:  "K5000": 

437:  if  B(K,4]=l;gto  "K5200" 

438:  fmt  "CALCULATE  NUMBER  OF  GUNS  REQUIRED" 

439:  I(K,2]/S(1,11*QI1,1] 

440:  I (K,41/S[l,2]-Q[l,21 
441:  UK, 21/3(2, 11*Q[2, 11 
442:  I [K,41/S[2,21-Q(2,21 
443:  for  L=1  to  2 
444:  for  N=1  to  2 
445:  int(Q [L,N1 )*t0 
446:  Q[L,Nl-r0*rl 
447:  if  rl#0;r0+l*Q(L,Nl 
448:  next  N 
449:  next  L 

450:  fmt  "CALCULATE  TOTAL  TIME  FOR  REPLENISHMENT." 

*7095 
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451;  for  L=1  to  2 

452;  for  N=1  to  2 

453;  T+TIL,N)-UIL,N1-»X[L,NI 

454;  next  N 

455;  next  L 

456;  frat  "CALCULATE  RATE  OF  FIRE." 

457;  for  N=1  to  2 
453;  X(1,N1*V(N,H1+1*R[1,N1 
459:  X(2,N1*W(N,H1+L>R(2,N1 
460:  next  t'j 

461:  fmt  "CALCULATE  TOTAL  NUMBER  OF  ROUNDS  REQUIRED." 
462:  for  L=1  to  2 


463; 
464: 
465: 
466 : 
467; 
468; 
469; 
470; 
471; 
472; 
473  ; 
47  4; 
475; 
476; 
477; 
478: 
479: 
480: 
481; 
482; 
483; 
434  ; 
435; 
486: 
437: 
488: 
43  9: 
490; 
4 91: 
492: 
493; 
494: 
495: 
496: 
497: 
493: 
499: 


for  N=1  to  2 

QIL,N1*R[L,N1-P[L,N1 

int(P[L,N]  )^r2 

P[L,Nl-r2-r3 

if  r3#0;r2+l-P[L,N] 

next  N 

next  L 

if  E$IL,3J="YES";gto  "K5100" 
fmt  4x,cl7;wrt  701,C$[I,J] 
fmt  " ";wrt  701 
"K5100"; 

1+17*1  ;J+17-»J 

fmt  " VOLUME  OF  FIRE  - HC  SMOKESCREEN" 

wrt  701 

fmt  " SCREEN  LENGTH  - METERS  = ",f3.2;wrt  701, R 

fmt  " SCREEN  DURATION  - MINUTES  = ",f3.2;wrt  701, T 

fmt  " ";wrt  701 

fmt  " GUNS  ROUNDS/MIN  TOTAL  ROUNDS" 

wrt  701 

fmt  " 105  ",f2.0,12x,f4.0,12x,f5.0 

wrt  701,0(1,11 ,V(1,H1 ,P[1,11 

fmt  " 155  ",f2.0,12x,f4.0,12x,f5.0 

wrt  701,Q[2,11 ,4(1, HI ,P[2,1I 
fmt  /;wrt  701 

fmt  •’  VOLUME  OF  FIRE  - WP  SMOKESCREEEN" 

wrt  701 

fmt  " SCREEN  LENGTH  - METERS  » ",f8.2;wrt  701 ,R 

fmt  " SCREEN  DURATION  - MINUTES  = ",f8.2;wrt  701, T 

fmt  " ";wrt  701 

fmt  " GUNS  ROUNDS/MIN  TOTAL  ROUNDS" 

wrt  701 

fmt  " 105  ",f2.0,12x,f4.0,12x,f5.0 

wrt  701,0(1,21 ,V(2,H1 ,P(1,21 

fmt  " 155  ",f2.0,12x,f4.0,12x,f5.0 

wrt  701,0(2,21 ,W(2,H1 ,P(2,21 
gto  "K5300" 

"K5200"; 


500;  1+16*1 ;J+16*J 
*20059 
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501:  fmt  " SMOKE  l^OT  REOOIRED  DUS  TO  ATMOSPHERIC  CONDITIOi'JS. " | 

502:  "K5300":fmt  /;wrt  701  C 

503:  next  K 

504:  "K5400'':  ^ 

505:  fmt  /;wrt  701 

506:  dsp  "DONE"  ’ 

507:  end  -> 

508:  "PNA":ret  exp(-r 6*pl/2) 

509:  "FNa":ret  axp(-pl“2)  | 

510:  "FNC":ret  exp(pi*r6*ln(.l/FlIl ) ) . * 

511;  "FND":ret  exp(-pl*r 6/4 . 1)  • 
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SAMPLE  CALCULATION 


KWIK  SMOKE  PROGRAM 


MET  SITE: 


ID 

DMN 

! LATITUDE  - DEG  = 

32.25 

1 LONGITUDE  - DEG  = 

107.72 

ALTITUDE  - KM 

1.32 

JULIAN  DATE  - DAY  = 

322 

ZULU  TIME  - HOUR  = 

23 

' METEOROLOGICAL  INPUTS: 

CEILING 

- METERS 

914.40 

CLOUD  COVER 

- PERCENT 

= 

40.00 

VISIBILITY 

- KILOIETERS 

= 

8.0  5 

PRECIPITATION 

3S 

NO 

TE.MPERATU  RE 

— DEG  C 

s 

10.00 

DEWPOINT 

- DEG  C 

= 

8.33 

WIND  DIRECTION 

- DEG 

= 

180.00 

WIND  SPEED 

- KNOTS 

= 

15.00 

AVE  ROUGHNESS  ELEMENT 

- CM 

= 

23.00 

METEOROLOGICAL  CALCULATIONS; 

PA30UILL  STABILITY  CATEGORY  = 
RELATIVE  HUMIDITY 


D 

89.37 


MONITIJN  EXPEUDITURES : 


VOLU.IE  OF  FIRE  - HC  SMOKESCREEN 
SCREEN  LENGTH  - METERS  = 

SCREEN  DURATION  - MINUTES  = 

200.00 

10.00 

GUNS  ROUNDS/ UN 

10  5 1 2 

155  1 1 

TOTAL 

ROUNDS 

17 

8 

VOLU  IE  OF  FIRE  - WP  SMOKESCREEEN 
SCREEN  LENGTH  - METERS 

SCREEN  DURATION  - MINUTES  * 

200.00 

10.00 

GUNS  ROUNDS/ 1 IN 

105  1 4 

155  1 2 

TOTAL 

ROUNDS 

43 

22 
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